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THE PROPAGATION OF AN ELECTROMAGNETIC WAVE ALONG 
AN INFINITE CORRUGATED SURFACE! 


By R. A. Hurp 


ABSTRACT 


Propagation of an electromagnetic wave along an infinite corrugated surface 
whose slot walls are vanishingly thin is investigated using a method based on the 
calculus of residues. The method yields expressions for the amplitudes and phase 
velocities which are highly accurate near cutoff when the number of slots per 
wavelength is greater than about five. 


1, INTRODUCTION 


In recent years, the application of corrugated surfaces to linear accelerators 
and antennas has aroused interest in the theoretical problem of finding the 
phase velocity and ‘“‘mode’”’ amplitudes of the surface wave which propagates 
on an infinitely wide surface. No rigorous solution to this problem is known, 
but approximate solutions have been obtained (1, 2, 4, 5). Of these papers, 
only that of Brillouin (1) has given expressions for the amplitudes, the others 
being concerned with the phase velocity only. 

The problem, as considered by Brillouin, was simplified by making the 
following assumptions: 

(a) The slot walls are vanishingly thin. 

(b) Only the propagating (TEM) mode in the slots is reflected from the 
base of the slots. As far as the higher order attenuated waves are concerned, 
the slots are of infinite depth.* 

To this idealized problem, Brillouin has given an approximate solution. 

In the work which follows, we shall be able to give an exact solution to this 
simplified problem, thus obtaining a useful solution to the problem of the 
corrugated surface with vanishingly thin slot walls, valid for frequencies 
near cutoff, when the number of slots per wavelength is about five or greater. 


2. PRELIMINARY THEORY 
Consider the configuration shown in Fig. 1. A perfectly conducting plate 
_ coincides with the plane y = —A. On the plate is built up a series of parallel, 
infinitely thin walls of height h, spaced a distance d apart. The walls are 


1Manuscript received July 9, 1954. 
Contribution from the Radio and Electrical Engineering Division, National Research Council, 


Ottawa, Canada, Issued as N.R.C. No. 3436. 
*Assumption (b) is equivalent to assuming that the slot width is small compared with both the 


slot depth and the free space wavelength. 
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Fic. 1. An infinite corrugated surface. 


assumed to extend to infinity in the +x directions, and to be infinite in number 
in the +2 directions. 

In the region above the slots it is assumed that there exists a wave proceeding 
in the +2 direction. This wave will be assumed to have three field components 
only, E,, E,, H,. The problem is thus a scalar one, with 


[1a] E, = —(it/we)dH,/02, 
[15] E,=  (t/we)dH,/dy, 
[2] ?H,/dy+H,/d2+hH, = 0. 


Here k = 27/) is the free space propagation constant. The time dependence of 
e“! has been omitted.* 

On account of the periodic nature of the structure, the field above the 
slots will not be a simple wave but will be a summation of an infinite number 
of waves of differing amplitudes and phase constants: 


[3] H, rea . ihe gente. 
[4] a, = (Bn? —k?)?. 


The phase constants, 8,, may be shown (1) to satisfy a relation of the form 
[5] Bn = Bot2an/d (n = 0, +1, +2,...). 


Since we are looking for surface waves, we shall assume that 6») > Rk and that 
the positive root is always taken in Equation [4]. 

In the vth slot region the field is represented by a superposition of all the 
x-independent modes which satisfy the boundary conditions at the base and 
sides of the slot: 

[6] H;” = DY Bn‘cos ¥m(y+h)cos(mxz,/d), 


[7] Ym = [R°— (ma/d)’}’. 


Here z, is measured from the left-hand wall of the vth slot, and is related to 
z through 


Z, = s—(v—3)d. 


*The m.k.s. unit system is employed. 
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The coefficients, B,“”, have to be chosen to give the correct phase change from 
slot to slot, hence 


B,, oad Be. 


The unknowns of the problem are the coefficient ratios A,/Bo, Bm/Bo, and 
the phase constant 8». These are determined by ensuring continuity of the 
tangential magnetic and electric fields (H, and 0H,/dy) across the mouth of 
the vth slot. Therefore 


[8a] : B,,cos(mrz,/d) cos ¥mh = > ha Cr. 
[8d] D YmBncos(mnz,/d) sin ymh = >> anAge #24, 
m=0 n=—c 


for 0 < z, < d. Solving Equations [8] for the coefficients Bm yields the equa- 
tions, true for all m, 





1 €mBmCOS Ymh _ 1BrAn genes 
[9a] 9 (— yr e -iBod alt x Qn a ’ 
d En¥mBn sin Ymh _ ss 1nByAne**!* 
[95] 9 (— \"e tee 1 — int Qn Fee ’ 
be s 
where En = J2, m 04 ‘ 
u, m>0o! 


Multiplication of Equation [9a] by tym and addition and subtraction of 
[9a] and [98] gives 


i8nd/2 





d €n¥mBme_* BnA ne” _ 9 
[10] 2 (—)Me ma ~ pp (m = 0,1,2,...), 
d cme BrAne oe = 9 
[11] 9 (— *, a—- L- an biy Ym (m == (), Re inate 


If now assumption () of the preceding section is invoked, it is seen that 
the factors e~™™* (m > 0) are small. We may therefore agree to neglect the 
left-hand side of Equation [10] when m > 0. Thus 


kdBo e —ikh > BrAne ore”? 


9 See emesis dil 
[12a] ge teed 1  “«,—ik ’ 
18,d/2 
[120] 0 = > BnA a pe oe 2% 6 ds 


Equations [12] are to be solved for the unknowns A». Equations [11] will be 
used later to determine the coefficients By. 
3. SOLUTION OF THE EQUATIONS 
Equations [12] have the same form as a set of equations solved by White- 
head (8), and the analysis used here follows his method. Consider a meromor- 
phic function f(w) which 
(i) has simple poles at w = a, (n = 0, +1, +2,...); 


| 
{ 
: 








af" 
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(ii) has simple zeros at w = tym (m = 1, 2, 3,...); 

(iii) tends to zero as |w| — ~, except near the poles an. 
Choose a circle C, which encloses the point w = 7k and the poles at w = an, 
|x| < s. Then by Cauchy’s theorem of residues 


uy lf Oe Fe poe 


830 201 Io, W—1Ym —ao An — VY; 


where r(a@,) is the residue of f(w) at w = an, and 
he 2 form = of 
*" ~ (0 for m > 0)* 


kdB — ikh 
[14] fk) = Sar, 


Setting 


the solution to Equations [12] is given by 
[15] A, = 1 (cn) Bn en”, 

A function of f(w) that satisfies all the conditions of the problem (i.e., 
conditions (i), (ii), (iii), and Equation [14]) is 


_ _AdBoe™™ ik—-ay Th (w)T2(ik) ‘ 
[16] fw) = — Sa] w—ay 1, (ik) Io) © exp [(d/ 7) (ik—w) log 2], 


where II,(w) and II;(w) are defined by 


I,(w) = [] (w—iyy)(—d/pr)e”™, 
p=1 


I(w) = [] (w—ay)(w—a_y) (d/2pm)*e”. 
p=l 


In Equation [16] the exponential factor is necessary to ensure the proper 
behavior of f(w) as |w| — ©. For details of the choice of this function, the 
reader is referred to Whitehead’s paper (8). 
4. CALCULATION OF THE MODE AMPLITUDES 

To compute the coefficients A», we should ordinarily make use of relations 
[15] and [16]. However, computation of the residues of f(w) as given by 
Equation [16] is awkward for the higher values of 7; so it will be found con- 
venient to derive an alternative expression. To do this, note that 


sin[d(w*+k*)*] 











[17] II, (w)I,(—w) = d(w*-+h?)' 
and 
d/2 ), Po— \W s da/2 0 hy 
[18] Ma(w)(—w) = SL /2) BoC ae sel e2 2) {Bot (w*+k*) 3] 
Hence 
4) = 1 kd'Boe“™” (ik—a0) (w+a0) Me(ik) T2(—w) 
[19] fw) = 4 é ipod aa (w° +k? i TI, (ik) TI, (—w) 
anes sin[d (w"- a ne & oo 2| 
» sin[(d/2){ Bo+ (w'-+k°)*} Jsin[(d/2) { Bo— (w' +R) }] (tk—w) log 2]. 
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Using Equations [15] and [19], we obtain, after some computation, 








A,| _ 1 k(antao) W2(—an) (ston) 
0) Bol ~ 2 alBl O(a) “PX? 


where use has been made of the relations 
|, (¢k)| = 1, 


sin (Bod /2) 
Bod /2 


which are easily derivable from Equations [17] and [18]. 

There does not seem to be any easy way of computing the coefficients 
Tl:(—a@,) and II.(—a@n) exactly. However, a good approximation is obtained 
by extracting the initial term from the infinite product and using the approxi- 
mations 


|Il2(¢k)| = 


a & |B,|, iy, = pr/d, 


in the remaining terms of the infinite products. The products can then be 
expressed in terms of gamma functions. Thus 








An| kd (antao) (@n+tar) (an-+a_1) 
aa Bo ~ Sr On| Brn| On +171 
x T (2+a,d/ 7) as | 24 log 3 
T[2+ (4/2) (an-+ Bo) P[2+ (d/2) (n— Bo] . - 


The coefficients |Ao/Bo|, |A1/Ao|, and |A_:/Ao| are shown plotted against 
h/X, for \/d = 20, in Fig. 2. 
For large n, 


[22] |An/Bo] = (kd/4)(m|n|)~°”. 





= eo 
o -05 10 15 
h/d 


Fic. 2. Variation of Ao/Bo, A-1/Ao, A1/Ao with h/d. 
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To compute the coefficients B», consider the expression 


pa SA ge a Nat 
[23] - 2 mt Ca wtiym la Beata +H Ym) 0, 
(m = 0, 1, 2,...). Comparing Equations [11] and [23], it is found that 
d ™ * 2 tymh 
[24] f(-t1m) = 5 Was nat 


Equations [16] and [24] then completely determine the Bn. For even m, 














, B, | _ _4ksin(Bod/2)_ |Ui(—iym)| E (: || 

: — = m\ — log 2—1 : 
aie) | Bo ~ dém|Ym| |i¥m-te| |IT2(—i¥m) | —" T = , 
and for odd m, 

‘i B,. 4k\cos(Bod/2)| |i (—tym)| i e )| 
9 = = sibs = - 
pr Bo dém|Ym| [?¥m-avo| [12(—7¥m) | — T log — 











When is large, 
[26a] |Bm/Bo| ~~ kd(xm/2)~*’sin (Bod /2) e*™", for even m; 


[26] |B,/Bo| ~ kd(am/2)~*”|cos(Bod/2) |e", for odd m. 


5. EVALUATION OF THE PHASE CONSTANT 


The phase constant, 8), may be determined by the elimination of By between 
Equations [14] and [24]. This leads to the relation 


tk—ao Wh(—ik)Mo(ik) 10] 2d og 2 | = nin 








ik+ao Th (tk) IIe (—ik) © 
This equation is satisfied identically in modulus for real 8). Upon taking 
arguments of each oi 


a onr 


45 2 eal 


Approximating the inverse sines by their arguments, and leaving out terms of 
order (d/\)* and higher, we obtain 


a if 94(1488) 4-00 ot(84) 420-2 al ca -1( 2), Hélon2 
me Bod Bo 


Tv 


[27] 





where ¥/(x) is the digamma function (3), 
= OY ss Bit L) 

and C = 0.577 is Euler’s constant. In isibcies: Eauntion [29] use has been 
made of the relation 
[30] V(x) +y¥(—x) = 2¥(1+x)+ mwcot rx —(1/x). 

Equation [29] can now be solved graphically for 8). This has been done in 
Fig. 3 which plots \,/A (Ag = 22/80) as a function of slot depth for various 
values of slot width. 
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Fic. 3. \,/A as a function of h/A for various values of \/d. 


6. DISCUSSION 


The results of the preceding sections call for some comment. With reference 
to Equation [29] it is apparent that for any hk which corresponds to a solution 
A,/X as given in Fig. 3 (say ho), there is an infinite set of values ho +pd/2 
(p = 0, 1, 2, 3,...) which corresponds to the same value of A,/A. The 
maximum value of / attained on any particular curve (Fig. 3) is the cutoff 
depth h,. For values of # greater than h, and less than \/2, it can be shown 
that it is necessary to take the a, megative, and hence there is no surface wave 
in this case. Summarizing, we have propagation when h is such that 


[31] (pA/2) <h < (pd/2)+h. (p = 0,1,2,...), 


and no propagation for / lying outside these limits. 

By differentiation of the right-hand side of Equation [28] it can be shown 
that the cutoff value of d/X, is 0.5. Then, from Equation [5], Bm = —B—m—13 
and from Equation [20], Am = A—m-—1. Thus a standing wave is set up on the 
surface. That this happens at cutoff has been pointed out by Walkinshaw 
(7). 

By substituting d/A, = 0.5 in Equation [29] and using the relation 

¥(3/2) = 2—C—2 log 2, 


one obtains 
(32] h, = (A/4)—(d/z)log 2 + (terms of order d?/)?). 


It is seen that, at cutoff, the slots may be considered to be parallel resonant 
with a discontinuity capacity at the mouth of the slots. The amount of short- 
ening of h, over \/4 given by Equation [32] is in agreement with that which 
would be expected from electrostatic considerations. (See Smythe (6).) 
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It is to be noted that the left-hand side of Equation [27] is invariant under 
the substitution of 1/d—1/\, for 1/A,. Hence Equation [27] has two roots 
A,’ and A,” connected by 1/A,’+1/A,” = 1/d. This is evident from the graphs 
in Fig. 3. Using Equation [5], Bm’ = —@—m-1’’; and it is apparent that the two 
solutions correspond to equal waves which proceed in opposite directions. 

When one compares the above results with those of Brillouin (1), it is seen 
that the expressions for the coefficients B», are considerably different. Thus 
we find 


[33] [Bn| = Ofm*"e-%™") for large m; 
whereas, according to Brillouin (in our notation), 

—tymh 

e 
[34] IBa| = lonfm—aa , for large m. 


It is not difficult to see that Equation [34] must be incorrect, since coefficients 
defined in this way lead to expressions for E, and E, which are finite 
everywhere. This is a violation of the edge condition which states that 
E,, E, = O{R™} for R—> 0, where R is the distance from the field point to 
any of the edges y = 0, 2 = (v+4)d. On the other hand it is relatively easy 
to show that our coefficients A,, Bm lead to expressions for E,, E, which 
exhibit the correct singularity at the edges.* 

It is possible to obtain better approximations to the solution of the general 
problem by keeping terms of order greater than zero in the left-hand side of 
Equation [10]. In order to establish the range of d/X and h/X for which the 
above results are valid, the calculation was carried out retaining the B, term. 
For \/d = 5 and 0.11 <h < h,, it was found that the new phase velocity 
differed by less than 1% from the value found above, while the amplitudes 
varied by less than 3.5°%. When \/d > 5, there is, of course, a correspondingly 
larger range of h in which comparable accuracy is obtained. 
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THE ROTATIONAL RAMAN SPECTRUM OF ETHYLENE! 


By J. ROoMANKO,? T. FELDMAN, E. J. STANsBURY,®? AND A. MCKELLAR‘ 


ABSTRACT 


The rotational Raman spectrum of ethylene excited by Hg4358 was photo- 
graphed with a reciprocal linear dispersion of 10.5 cm.~! per mm. In addition to S 
and R branches which, in the symmetric top approximation, can be used to deter- 
mine Bfo} + Cfo}, the spectrum shows a series of lines which must be interpreted 
in terms of the asymmetric top model. This series is shown to be due to the tran- 
sitions, AJ = 2, between the two lowest 7 levels of each rotational state, and is 
used to determine Co}. The rotational constants obtained from the analysis are 
Ajo} = 4.66+0.20, Bio) = 1.008;+0.006, Cio) = 0.828 5+0.003 cm. 


INTRODUCTION 


The plane symmetric form (point group V,) for the ethylene molecule is 
generally accepted, and the structures of the vibrational infrared bands and 
the rotational Raman band show that the molecule is very nearly a symmetric 
top. If the principal moments of inertia are designated 74, Jz, and I¢, with 
I, < Ig < Ic,* then Ig refers to the principal axis of inertia perpendicular 
to the plane of the molecule, and J, and J, to the axes in the plane, which are 
parallel and perpendicular, respectively, to the carbon-carbon bond. Since 
I, is nearly equal to Jg the symmetric top approximation has been used to 
determine the values of the rotational constants A, B, and C (A = h/8r°cI4, 
etc.) from the spectrum. The value of B = 3(B + C) can be obtained from 
the P and R branches of parallel infrared bands and the value of A — B 
from the Q branches of perpendicular bands. Using the relationship for planar 
molecules 1/C = 1/A + 1/B, the constants A, B, and C and the moments of 
inertia can be calculated. The most accurate determination of these constants 
from the infrared spectrum appears to be that of Gallaway and Barker (3), 
who obtained B;,; = 0.9116 and A;,; = 4.867 cm.— for the ground vibrational 
state. 

The constants Aj}, Bro, and C,.; can be also determined from rotational 
spectra. However, since the molecule has a center of symmetry, a pure rota- 
tional spectrum exists only in Raman effect. The Raman spectrum was 
photographed by Lewis and Houston (6) using a Hilger El quartz spectro- 
graph and Hg2537 as the exciting radiation. With the available dispersion 
(~43 cm.~! per mm.) only a single series of lines with an average spacing of 
3.68 cm.—! was observed. The series was interpreted as the S branch (AJ = 2) 
in the symmetric top approximation and the value B;,; = 0.920 cm.~! was 
calculated from the average spacing. This value is in fair agreement with the 
infrared result. 


1 Manuscript received August 30, 1954. 

Contribution from the McLennan Laboratory, University of Toronto, Toronto, Ontario, and to 
be reprinted as Contributions from the Dominion Astrophysical Observatory No. 40. 

2Holder of a Garnet W. McKee - Lachlan Gilchrist Scholarship, School of Graduate Studies, 
University of Toronto, 1952-54. 

3Present address: Bell Telephone Laboratories, Murray Hill, New Jersey, U.S.A. 

‘Dominion Astrophysical Observatory, Victoria, B.C. Visiting Professor of Astronomy and 
Physics, University of Toronto, 1952-68. 

*The nomenclature used follows that of Herzberg (4). 
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The development of sources of high intensity (11, 9, 12) has made it possible 
to photograph Raman spectra with much higher dispersion than was formerly 
possible. In the present investigation the rotational spectrum of ethylene has 
been re-examined with the aim of improving the accuracy of the molecular 
constants and of ascertaining whether departures from the symmetric top 
spectrum might be observed in high dispersion spectrograms. 

EXPERIMENTAL 

The multiple reflection Raman tube and the high dispersion spectrograph 
have been described (12). The experiments with ethylene were carried out in 
an early stage of the development of the equipment when the irradiating 
source was one mercury lamp of the shape described in Reference (10), with 
internal cooling in the discharge section. Research grade ethylene gas was 
used at pressures from 1.5 to 3 atm. Exposure times from 2 to 24 hr. were used 
to record the spectrum on Kodak 103a-0 spectroscopic plates with a spectral 
slit width of 0.33 cm~!. The reciprocal linear dispersion of the spectrograph 
was 10.5 cm.~! per mm. The microphotometer traces reproduced in Fig. 1 
show the Stokes portion of the spectrum for exposure strengths of 36 and 22 
atm. hr. with Hg4358 as the exciting line. 


S' Branch ie ee ee 
: (a) 
o 
wy 
ao 
vr 
wo 
oh 
R Branch (b) 
J?6 ' iO] [15] 
jc SS SO CIS oon) Se gee ee man 
S Branch 
Got VP ee 0° 75 100. 


Raman Displacement (cm.") 


Fic. 1. Microphotometer traces of the Stokes portion of the rotational Raman spectrum of 
ethylene with exposure strengths of (@) 36 and (6) 22 atm. hr. 


An iron arc comparison spectrum was photographed in juxtaposition with 
each Raman spectrum. Since there is a scarcity of iron lines in the region of 
Hg4358 a lightly exposed rotational Raman spectrum of oxygen was also used 
for comparison. The frequency shifts of the oxygen lines were calculated from 
the very accurate molecular constants given by Babcock and Herzberg (1). 
The wavelength of the exciting line was taken as 4358.35 A as given in the 
M.I.T. Wavelength Tables. Three spectrograms with different strengths of 
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exposure, 36, 28, and 5 atm. hr., were measured on a Hilger comparator. The 
Stokes and anti-Stokes frequency shifts were averaged for each plate; finally, 
the frequency shifts for the three plates were averaged. The results are given 
in Table I. 























TABLE I 
ROTATIONAL RAMAN SHIFTS OF C2H, 

Av (cm.~') _ Intensity* Assignment Av (cm.~') _ Intensity* Assignment 
7.27 5 R(3) 59.32 6 S’(16) 
9.00 5 S(1), R(4) 60.42 7 S(15) 

10.60 3 — 62.73 2 S’(17) 
11.67 3 — 64.39 +4 S(16) 
12.85 7 S(2), R(6) 66.27 4 S’(18) 
14.64 3s R(7) 67.95 3 S(17) 
16.49 8 S(3), R(8) 69.88 1 S’(19) 
18.27 4 R(Q) 71.78 3 S(18) 
20.15 9 S(4), R10) | 72.79 3d S’ (20) 
21.99 3 R(11) 73.79 1 — 
23.80 8 S(5), R(12) | 75.56 3 S(19) 
25.70 3s R(13) 79.32 3 S(20) 
27.54 8d S(6), R(14) 80.37 Is S’ (22)? 
29.35 3s R(15) I] 82.83 3d S(21) 
31.15 8 S(7), R(16) 86.35 3d S(22) 
33.02 4 S’(8) 89.82 2dd S(23) 
34.99 a S(8) | 92.7 2 S’ (26) 
36.07 4 S’(9) 93.68 2 S(24) 
38.51 e S(9) i 96.95 ld 5(25) 
39.75 4 S’(10) | 99.40 ls S’(28) 
42.43 7 S(10) | 100.92 ld S(26) 
43.39 4 S’(11) | 102.56 ls 5S’ (29) 
46.10 7d S(11) 105.18 ldd (27) 
49.72 7 S(12) | 109.02 1d S(28) 
51.78 1 S’(13)? || 112.38 1d S(29) 
53.29 6 S(13) | 115.97 Od S(30) 
56.80 6 S(14) 118.78 Od S(31) 


*The intensities given are visually estimated on a scale 10 to 0. Sharp lines are designated by s, 
diffuse lines by d, and the more diffuse lines by dd. 


DISCUSSION OF THE RESULTS 

The S lines, corresponding to the transitions AJ = 2, AK = 0, in the sym- 
metric top approximation, are easily recognized in the spectrum. These were 
recorded by Lewis and Houston (6) up to J = 20 and, in the present investiga- 
tion, to J = 31. The R lines (AJ = 1, AK = 0) up to J = 15 are also easily 
observed on the high dispersion spectrograms. The relatively low intensity 
of the R branch is due to the low value, ~0.25, of the ratio J,4/I, (Placzek 
and Teller (8)). The assignments for the lines of the S and R branches are 
unambiguous and are listed in Table I. 

In addition to the R and S lines other lines are observed in the part of the 
spectrum where the R branch is very weak. These extra lines form a series 
with a spacing somewhat smaller than that of the S lines and will be referred 
to as the S’ lines. The S’ series, observable only when it is out of step with the 
S series, imparts to the spectrum the characteristic appearance shown in 
Fig. 1. 
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When the Raman displacements of the S lines are plotted against J, as in 
Fig. 2, a straight line can be drawn through all the points. Also, when the 
correct relative numbering of the lines has been determined, the displacements 
of the S’ lines plotted against J yield a straight line. The absolute numbering 
of the S’ series in Fig. 2 is based on the analysis given below. 
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Fic. 2. The S and S’ frequency shifts in the rotational Raman spectrum of ethylene as a 
function of the rotational quantum number J. 


The R and S branches of the spectrum can be interpreted in terms of the 
symmetric top approximation. The S’ branch, on the other hand, is an indica- 
tion of the deviation of the molecule from an exactly symmetric top, and its 
interpretation must therefore be based on the theory of the asymmetric top. 

In the symmetric top there are J + 1 sublevels of different energies for each 
value of the total rotational angular momentum quantum number, J; the 
assembly of energy levels is given by the expression 
[1] F(J, K) = BJ(J +1) + (A — B)K*, J=0,1,2,..., 
where K = 0, 1, 2,...J and each level with K # 0 is doubly degenerate. In 
the asymmetric top this degeneracy is removed; for each value of J there are 
2J + 1 sublevels, F(J-), which are ordered by the subscript 7 taking the values 
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—J,-J+1,...+J-—1,+J with r = —J assigned to the lowest level. 
When the asymmetric top is very nearly a symmetric top the 7 levels can be 
correlated with the K levels of the equivalent symmetric top. Thus, for the 
prolate case, Iz ~ Ic, the highest pair of 7 levels corresponds to K = J and 
the lowest 7 level to K = 0. 

To determine the origin of the S’ lines the frequencies for the J = 12 <— J = 10 
Raman transition were calculated using the rotational constants as determined 
by Gallaway and Barker (3) and the eigenvalues of the asymmetric top as 
evaluated by King, Hainer, and Cross (5). The populations of the various 
initial levels were calculated using the statistical weights as given by Wilson 
(13). Since the matrix elements of the polarizability have not been evaluated 
for the asymmetric vop, the matrix elements for the symmetric top (8) were 
used to determine the intensities. This procedure restricts the possible tran- 
sitions to those for which AK = 0 in the equivalent asymmetric top; however, 
since the ethylene molecule does not deviate much from the symmetric top 
model, the calculated spectrum should not be greatly in error. 

In Fig. 3(a) the spectrum of the 12-10 transition is plotted for the case of 
very high spectral resolution. In Fig. 3(6) the spectrum is modified for a line 
width of 0.5 cm.~!, which corresponds approximately to the resolution in the 
experiment. The highest intensity is that due to the nearly superposed lines 
resulting from transitions at high r-values; since these energy levels deviate 
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Fic. 3. (a) Calculated structure of the J = 12<—J = 10 transition in the rotational 


Raman spectrum of ethylene. 
(6) The calculated structure as modified by a line width of 0.5 cm™. 
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very little from those for the equivalent symmetric top this maximum corre- 
sponds to the S Jine. The next highest intensity is that given by the two transi- 
tions involving the two lowest sublevels in each J state (J’, =12_32- J", 
= 10-19 and J’,,=12_1:<- J”. =10_9); for these levels the eigenvalues 
deviate farthest from those of the symmetric top. It is therefore evident that 
the S’ lines can be correlated with the transitions, AJ = 2, between the two 
lowest sublevels of each J state. As has been shown by Dieke and Kistiakowsky 
(2), these two sublevels approach each other at high J values. The minor 
maxima in Fig. 3(6) are not observed, probably because of overlapping with 
lines from other transitions in J. 

For the sublevels with the highest values of 7 the formula for the equivalent 
symmetric top, 


[2] Foroate(J, K) = 3(B + C)J(J + 1) + [A — 3(B + C))K?, 


is an adequate representation of the energy. Applying the selection rules, 
AJ = 2, AK = 0, one obtains for the Raman shifts in the Stokes and anti- 
Stokes S branches the expression 


(3] [Ar] = 2°0B+C)J+3(B+C), J=0,1... 


The graph of |Av| against J is thus a straight line with slope 4B = 2(B + C). 
As is evident in Fig. 2, the experimental data for the S lines can be represented 
by a straight line over the whole range of J. Since there is no noticeable curva- 
ture in the graph it is certain that centrifugal stretching terms, not taken into 
account in [2], can be neglected. A least squares analysis of the data gives the 
value B,,; = 0.9187+0.0017 cm™. 

In the proposed interpretation of the S’ lines an analytical representation of 
the S’ frequency shifts requires an expression for the average value of the 
energies of the two lowest 7 sublevels. For a molecule which approaches the 
limiting case of an oblate symmetric top Mecke (7) obtained the approximate 
formula 


[4] 3[F(J-2)+F(J—s1)] = CI°+3(A+B)J 


aes _) ( .) 
1(A Bye Wits) 


B~A (=0.910 for C2H;). 


2(C—3(A+B)] 


where b* = 


Numerical calculations were carried out to ascertain the degree of accuracy to 
which this formula can be used for ethylene, which approaches the limiting 
case of a prolate symmetric top. In Table II the values of 4[F(J_7) + F(J—741)] 
are given for the even values of J up to J = 12. ‘“‘Correct”’ values calculated 
from the eigenvalues of King, Hainer, and Cross are compared with values 
calculated from formula [4]; in both cases the rotational constants as deter- 
mined in this investigation (Table IV) were used. The differences, 6, between 
the ‘correct’ and the approximate values are given in the last column in the 
table. It is apparent that, although 6 increases with J, the differences between 
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TABLE II 
CALCULATED VALUES OF }3[F(J_y) + F(J_741)] 





Approximate value 











SF “Correct’’ value calculated from [4] Difference, 6 
(cm.~) (cm.~) (cm.~) 
4 19.783 19.914 0.131 
6 39.245 40.219 0.974 
8 65.672 67.117 1.445 
10 98.780 100.633 1.853 
12 138.584 140.775 2.191 





successive 6’s decrease with J; this indicates that formula [4] might be used to 
calculate the frequency shifts of the S’ lines at high J values. 

With the selection rule AJ = 2 the frequency shifts of the S’ lines using 
formula [4] are 


[5] |Av| =" 3{F[(J+2)_5-2]+F[(J+2)_s-1]} —${F(J_s)+F(J_341)} 


= 4C(J+1)+(A+B)-3(4 ~sy(14+5)(4-2,) ,J=0,1,2.... 


In Table III ‘‘correct’’ frequency shifts for various transitions are compared 
with those calculated from the approximate formula [5]. The differences, A, 


TABLE III 


Approximate value 





Transition “Correct”’ value calculated from [5] Difference, A 
af (cm.~') (cm.~!) (cm.~) 
64 19.462 20.305 0.843 
8—6 26.427 26.898 0.471 

10<—8 33.108 33.516 0.408 
12—10 39.804 40.142 0.338 





decrease with increasing J, but the rate of decrease is less marked at high J 
values. However, for the 1210 transition the discrepancy is already of the 
order of 0.3 cm™. It is therefore concluded that Mecke’s approximation for 
the oblate asymmetric top can be safely applied to the case of ethylene. 

Since the S’ lines are observed only for J > 7 where they are no longer 
obscured by the R lines, a simplification in formula [5] can be made. With 
increasing J the term in 1/(J? — 1) decreases rapidly so that, at high J values, 
the S’ frequency shifts plotted against J should give a straight line with the 
slope 4C. A least squares analysis of the data, plotted in Fig. 2, gives the value 
Cio) = 0.8289+0.0027 cm™!. It may be noted that the straight line in the plot 
‘of |Av| vs. J at high J values has an intercept on the |A»| axis equal to 


4C + (A + B) — 3(A — B)d*(1 + B*?/4) 
according to formula [5]. This fact was used to obtain the absolute numbering 
of the S’ lines given in Table I and Fig. 1. 
The rotational Raman spectrum of ethylene at high dispersion can thus be 
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used to evaluate all three rotational constants of the ground vibrational state. 
From the experimental values of By.) = 1/2(Bto) + Cto)) and Ct, and the 
relation 1/A = 1/C — 1/B, for a planar molecule, the values of Aj.) and 
Bi.) can be obtained. The results with their probable errors are summarized 
and compared with the infrared values of Gallaway and Barker (8) in Table IV. 
TABLE IV 
ROTATIONAL CONSTANTS FOR THE GROUND VIBRATIONAL STATE OF C2H, 








From the From the 
Raman spectrum __ infrared spectrum (3) 
(cm.~) (cm.~) 
Bio} 0.9187 +0.0017 0.9116 
A [0] 4.66 +0.20 4.867 
Bio} 1.008; +0.006 0.996 
Cio} 0.828, +0.0027 0.827 


The value of C,,; found in the present investigation agrees very well with the 
infrared value. The agreement of the values of B,,) is not as good, but the 
Raman value is probably the more accurate. The Raman value of Aj; has a 
large probable error because it is calculated from the difference of two quan- 
tities which are almost equal in the relation 1/A = 1/C — 1/B. The infrared 
value of A;,}, determined from the spacing of Q branches in perpendicular 
infrared bands, is undoubtedly more accurate. 

With further advances in the techniques of high resolution Raman spectros- 
copy it is quite possible that more detail of the rotational Raman spectrum of 
ethylene can be resolved with a resultant improvement in the molecular 
constants. Also, the analysis of the rotational Raman spectra of slightly 
asymmetric top molecules, based on the observation of S and S’ series, may 
well find applications in heavier molecules. 
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STATISTICAL STUDIES OF POLAR RADIO BLACKOUTS! 


By J. W. Cox AND KENNETH DAVIES 


ABSTRACT 


A statistical study of high frequency radio blackouts in Canada is made from 
records taken at several ionosphere sounding stations. Both vertical incidence and 
communication data are examined to determine the geographical, seasonal, and 
diurnal distributions of the frequency of occurrence of blackout. It is found that 
blackouts are most abundant in the morning and that the time of maximum 
occurrence increases with increasing latitude. 


I. INTRODUCTION 

It is well known that high frequency radio communication is more difficult 
at high latitudes than at low latitudes because of the frequent occurrence of 
ionospheric disturbances. Studies by means of conventional vertical incidence 
pulse soundings have shown that, in general, the same ionospheric regions 
exist at both high and low latitudes but that at high latitudes rapid variations 
of considerable size in the ionospheric characteristics are more frequent and 
pronounced. This is particularly true of the F, and D regions which are the 
regions of greatest importance for high frequency radio propagation. 

In the present paper a description will be given of some statistical studies 
of the frequency of occurrence of the so-called ‘‘no-echo”’ condition at a number 
of stations in Canada and in the U.S.A. over the period 1949 to 1952. Few 
attempts appear to have been made to study blackout phenomena and the 
present study has been made in an effort to make some use of the available 


data. ’ 
II. SOURCES OF DATA 


The material used in this study consists of: the hourly records of observa- 
tions at appropriately chosen ionospheric sounding stations at which records 
of vertical incidence pulse soundings are made, and the reception, at thesé 
ionospheric sounding stations, of signals from WWYV, the standard high fre- 
quency sender of the National Bureau of Standards, Washington, U.S.A. 
These records were not made primarily for the purposes of this study, and 
certain limitations on accuracy, discussed later, had to be accepted. 

The locations of the various stations are shown in Fig. 1 and in Table I. 


TABLE I 
POSITIONS OF IONOSPHERE SOUNDING STATIONS 





? r & 

Station North West North L.M.T. 
Resolute Bay ae 94.9° 82.9° 90° 
Baker Lake 64.3° 96.0° 73.7° 90° 
Churchill 58.8° 94.2° 68 .8° 90° 
Winnipeg 49.9° 97.4° 59.6° 90° 
Washington 39.0° 77.5° 50° 75° 


@ Geographical latitude. 
X Geographical longitude. 
& Geomagnetic latitude. 
1Manuscript received July 12, 1954. 
Contribution from Defence Research Telecommunications Establishment, Defence Research 
Board, Ottawa. Research conducted under Defence Kesearch Board Project D48-95-1i-02. 
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Fic. 1. Locations of ionosphere sounding stations. 


(a) Vertical Incidence Results 

Since there are no quantitative measurements of ionospheric absorption 
available for the stations, attention has been confined only to those occasions 
when no echoes were received. These occasions are denoted by the symbol B 
on the data sheets, and in the remainder of this paper, the occurrence of 
blackout at vertical incidence is synonymous with the occurrence of ‘‘B’’. A 
total loss of echo means a B on all data sheets, and apart from occasional 
errors, all sheets should give identical results. In fact the fin sheets have been 
used for all stations except Washington for which, as the fain was not available, 
the foF: sheets were used. 

This approach is qualitative and does not show how much more power 
would be required to produce a visible echo, nor does it allow for possible 
influence of noise variations. The disturbances dealt with are known to be of 
considerable magnitude, and are essentially an on-off phenomenon. It is, 
therefore, felt that a considerable change in equipment performance, or con- 
siderable differences between stations, would not have a great influence on 
the results obtained. 

The observations as taken at the stations record only 15 sec. of each hour. 
This makes studies of blackout duration precarious, but provides enough 
material for a statistical study of the probability of finding a blackout in 
progress at any random moment of sampling. 

(b) Reception of WWV 

Subjective measurements of the strength of signals received from WWV 

were also made hourly. These measurements were of no use as attenuation 
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measurements and, therefore, an on-off criterion similar to that used in the 
vertical incidence studies was adopted. At all stations, at least one of the trans- 
missions on 2.5 Mc., 5.0 Mc., 10.0 Mc., and 15.0 Mc. was normally audible. 
A blackout is considered to have occurred when none of these frequencies 
could be heard at all. 


(c) Effect of foF2 Variations 

Before discussing the results in detail one fact should be emphasized which 
may distort our qualitative picture. It is to be expected that the attenuation is 
produced by ionization in the lower part of the ionosphere, and that the 
absorption will vary roughly inversely as the square of the frequency (see 
Appleton (1)). Thus, for a given amount of ionization, corresponding to a fixed 
absorption on any one frequency, there is a greater likelihood of recording a 
blackout at times of low foF2. The latter happens diurnally, seasonally, and 
with the long period solar activity variation. It will be shown that the observed 
frequency of blackouts does not change in sympathy with the f,F:2 variations, 
and it is, therefore, considered that though there may be a small contribution 
from this cause, the blackout phenomenon is of a sufficiently ‘‘on—off’’ nature 
that the results would not be greatly altered if foF2 variations were allowed for. 
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Fic. 2. Occurrence of blackouts as a function of latitude. Wn. Washington, Wg. Winnipeg, 
C. Churchill, B. Baker Lake, R. Resolute Bay. 
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III. ANALYSIS OF VERTICAL INCIDENCE DATA 


(1) The Geographical Distribution of Blackout 

In Fig. 2 the totals of ‘Bs’ recorded at the stations over the four-year period 
1949 to 1952 have been plotted against geographical latitudes. This indicates 
that the occurrence of blackout is greatest near Churchill. 
(2) Temporal Variations 

(a) Secular Variations 

The long term variations in the frequency of occurrence of blackout pheno- 
menon are shown in Figs. 3 and 4 which show the yearly totals and the indi- 
vidual monthly totals of hours of blackout respectively, independently of time 


of day. 
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Fic. 3. Secular variation in frequency of occurrence of blackout. 


Comparison of the blackout data with the Zurich relative sunspot number, 
also plotted in Fig. 3, shows that peak values of ‘B’ occurred in 1947 at 
Churchill, suggesting some dependence of blackout phenomena on solar 
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activity. However, the data are too scanty to draw a definite conclusion, 
especially as the data for Baker Lake and Resolute Bay appear to have no 
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Fic. 4. Monthly totals of B. 
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Fic. 5. Monthly mean totals of B (1949 to 1952). 
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connection with sunspots. A more interesting correlation is that between the 
occurrence of blackout and the mean annual values of the magnetic Kw figure, 
which are also plotted in Fig. 3. The correlations with blackouts at Churchill 
and Winnipeg are fair and the correlations with the data for Baker Lake and 
Resolute Bay give some support to the conclusion that, at least statistically, 
blackouts are associated with enhanced magnetic activity. The lack of detailed 
correlation may be more significant than the slight over-all correlation. 

(b) The Seasonal Variations 

The individual monthly totals of ‘B’ are plotted in Fig. 4. There is a slight 
tendency for peak values to occur near the equinox months at the two southern- 
most stations. This is brought out in Fig. 5 in which the mean monthly totals 
for the four years 1949, 1950, 1951, and 1952 are plotted. There is very little 
seasonal variation at Baker Lake while at Resolute Bay there is a well-defined 
increase from winter to summer. 

(c) The Diurnal Variations 

In Fig. 6 the total numbers of ‘Bs’ occurring at each hour of the day over the 
years 1949 to 1952 are plotted. These graphs indicate a well-defined diurnal 
variation. The chief points to be noted are: 


CHURCHILL 





Fic. 6. Diurnal variation of frequency of occurrence of B. 
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(i) Daily maxima of occurrence in the morning, except at Resolute Bay, and 
minima of occurrence in the evening. A further study of the Resolute Bay data 
(not reproduced here) shows that the early morning maximum at Resolute 
occurs almost entirely in the winter months. 

(ii) Although the total number of hours of occurrence of blackout at 
Washington is too small to construct a reliable diurnal curve, it is clear that 
most blackouts occur there during the hours of darkness. It may be inferred 
from this that the effects of sudden ionospheric disturbances of the Dellinger 
type are small since nearly the same number of Dellinger blackouts occur at all 
stations, except in the winter when the northern ones may be in darkness. 
Over the four-year period 1949 to 1952, 89 blackouts were recorded at Washing- 
ton; 66 of these occurred between 00:00 hours and 06:00 hours. Thus even if 
all the remaining 23 were Dellinger type fadeouts they would constitute a very 
small fraction of the total numbers of blackouts observed at the Canadian 
stations. 
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Fic. 7. Plot of times of diurnal maxima of occurrence of blackout against latitude. 


(iii) A plot of local mean time of maximum occurrence against latitude 
gives a remarkably linear relation (Fig. 7) to which the second maximum at 
Resolute Bay fits quite well. The relation is 

T = -—5.7+ 36 


when T is local mean time (hours) and ¢ is geographical latitude in degrees. 
This variation which is one hour in four degrees of latitude recalls the work of 
Stagg (3) on the diurnal variation of magnetic disturbances. He found a maxi- 
mum in the late evening but states that ‘‘there is a disposition for the time of 
the maximum to be delayed with increasing latitude at the rate of approxi- 
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mately one hour every 5 
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Fic. 8. Diurnal variation of the minimum frequency reflected from the ionosphere (hourly 
median values 1949 to 1952). 


(iv) In order to see whether the diurnal variation in the occurrence of ‘B’ is 
brought about by the normal diurnal variation of D-region absorption the 
hourly median values of the minimum frequency reflected from the ionosphere 
have been plotted in Fig. 8. Now, provided the equipment characteristics are 
constant, finn is largely determined by ionospheric absorption. From Fig. 8 
it will be seen that the maxima are attained near noon at all the stations. It 
may be concluded, therefore, that the observed diurnal variation is not 
produced simply by the variation in the position of the sun. It should be noted 
that the observed diurnal variation of fain is diminished because the receiver 
gain is greater by day than by night but this would not affect the time of 
maximum value of fin. 


(3) Duration of Blackout 


It is of some interest to examine the distribution of the total time of blackout 
amongst the individual blackouts of various durations. To study this aspect 
of the problem the data have been classified according to the number of 
successive hourly records to which the symbol ‘B’ is assigned. In Fig. 9 the 
distribution of the total blackout time has been plotted against the duration 
of the blackout. It is inherently assumed that if two successive hours are 
blacked-out the whole of the interval between is also blacked-out. Hence, the 
occurrence of the symbol ‘B’ on n successive hours is interpreted as representing 
a blackout of length longer than m — 1 hours but shorter than m + 1 hours. It 
may be relevant to note here that a study of short duration blackouts suggests 
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Fic. 9. Distribution of the total time of blackout amongst blackouts of different durations. 


that the influence of sudden ionospheric disturbances of the Dellinger type is 
relatively slight. 

From Fig. 9 it can be seen that the distributions of the total time of blackout 
amongst the individual blackouts are very similar at Winnipeg, Churchill, 
and Resolute Bay. There is a remarkable similarity in the relative distributions 
of short duration blackouts. Thus it appears that although a blackout may be 
less likely at Resolute Bay than at Churchill, the probability that it will last 
up to about 24 hr., say, is roughly independent of the locality. It appears, 
however, that there is a much larger percentage of one-hour blackouts at 
Baker Lake than elsewhere. 


(4) Extent of Blackouts 

It is seen that many more blackouts occur at some stations than at others, 
suggesting that such blackout phenomena are to a large extent localized in 
space. It is of considerable interest to study to what extent the occurrence of 
blackout at one station is reproduced at another. It is not possible in a reason- 
able space to consider all possible correlations and so attention will be confined 
to a detailed correlation between blackouts which occurred at Churchill and at 
Baker Lake for 1951. 

Blackouts of all durations will be considered. For the present purpose it is 
convenient to express the correlation as the percentage of Baker Lake black- 
outs which occurred simultaneously at Churchill. The monthly values thus 
obtained are plotted in Fig. 10. In considering such correlation it must be 
remembered that even if the occurrence of blackouts were completely random 
at both sites, appreciable correlation would still be evident. Values to be 
expected in the case of random distributions can be estimated from the known 
frequencies and durations at the two stations, and are also included in Fig. 10. 
From this diagram it is clear that the observed values are significantly greater 
than would be expected on the basis of chance. 
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It should be emphasized that most of the simultaneous occurrences of 
blackouts are due to long duration disturbances. Thus a study of blackouts of 
duration of one hour at all four stations for 1951 showed that there was no 
occasion on which a blackout of one-hour duration occurred at all four stations 
simultaneously. 

It thus appears as if blackouts of short duration are essentially localized 
phenomena and that those of long duration are more widespread. 


IV. RECEPTION OF WWV 
(1) Geographical Distribution 
To illustrate the geographical distribution the total number of hours of 
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blackout have been plotted against geographical latitude for the four Canadian 
stations. This is shown in Fig. 11. The variation of these blackouts with 
latitude bears a striking resemblance to that of the vertical incidence case 
shown in Fig. 2. There are two possible explanations for this similarity, namely, 
(i) that reception of radio signals is governed by the state of the ionosphere 
near the receiving site, 
or (ii) that both vertical incidence soundings and reception of long distance 
signals depend on operational conditions at the stations, e.g., noise level, 
operating procedure, etc. 
The former explanation is very unlikely, since normally it is assumed that 
reception of long distance signals is governed by conditions along the propaga- 
tion path. On the other hand the alternative explanation is quite likely, 
because both types of observations depend to some extent on the personal 
judgments of the observers. 


(2) Temporal Variations 

(a) Secular Variations 

In Fig. 12 are plotted the total number of hours of blackout independent of 
time of day. It is clear that the long term variations at Winnipeg and at 
Churchill are in an opposite sense to the sunspot cycle variation. It appears, 
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Fic. 12. Number of hours of no reception of WWV per month on frequencies of 2.5, 5.0, 
10.0, 15.0 Mc. 


therefore, as if these blackouts are due to electron limitation effects in the E 
and F regions. The data for Baker Lake and Resolute Bay are too meager to 
reach any definite conclusions but the frequency of occurrence of blackout at 
these stations appears to vary in phase with the sunspot cycle. 


(b) Seasonal Variations 
The mean monthly totals of hours of blackout for the years 1949 to 1952 
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Fic. 13. Mean monthly totals of hours when WWV was not received 1949 to 1952. 


inclusive are plotted in Fig. 13. While there appears to be little definite seasonal 
variation, it may be relevant to note the very low July values which occur at 
all four stations. 


(c) Diurnal Variations 

The diurnal variations of the frequency of occurrence of blackout are shown 
in Fig. 14 where the hourly totals over the four-year period 1949 to 1952 are 
plotted. Here again there is a marked similarity with the vertical incidence 
phenomenon, because most blackouts occur before noon, and the ‘‘double 
hump” is again evident at Resolute Bay. It is interesting to note, however, 
that the times of maximum frequency of occurrence are between one and two 
hours earlier than in the case of the vertical incidence curves. This is to be 
expected qualitatively from the geometry of the paths and from Fig. 7. 


(3) Duration of WWV Blackouts 

The data have been classified in the same way as for the vertical incidence 
data—see Section III (3). The distribution of the total time of blackout 
amongst the individual blackouts is shown in Fig. 15. The results are similar 
to the vertical incidence results. The probability of occurrence of short period 
blackouts is about the same at all four stations. Long period blackouts appear 
to be more frequent in the north than in the south. 
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V. SUMMARY 


Although the data used for the present studies are somewhat uncertain 
because of their subjective nature, it is believed that the results will be of 
interest in a general way. Owing to the difficulty of making any accurate 
ionospheric absorption measurements in the Arctic regions, very little work 
has been done on radio wave attenuation at high latitudes. The present 
studies, therefore, constitute an attempt to obtain some useful information 
from available data. 

Summarizing the most important results: 

(i) For about 5 to 10% of the total time no echoes are obtained from the 
ionosphere at vertical incidence, 

(ii) Blackouts usually occur in the morning hours, 

(iii) The time of daily maximum of occurrence of blackouts is roughly a 

linear function of latitude, 

(iv) There is no well-defined seasonal variation common to all stations, 

(v) Short duration blackouts appear to be localized in space. 

(ii), (iii), and (v) suggest that blackouts are caused by corpuscular action 
(Piggott (2)), rather than by ultraviolet radiation. 

Probably the most interesting results are (ii) and (iii). The cause of the 
morning maximum and its dependence on latitude is not known. It is of course 
obvious that there must be some observer contribution to this essentially 
qualitative work, and it may be that it is significant. It is for this reason that 
the details of the apparent variation of blackout occurrence with latitude are 
not greatly emphasized. It is not believed that observer judgment would 
produce all the features obtained from the records and therefore these studies 
are presented to show that there certainly are phenomena worthy of more 
detailed and deliberate study. 
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SOME MASS SPECTROMETRIC DATA ON PHOSPHORUS! 


By LARKIN KERWIN 


ABSTRACT 


Investigation of red phosphorus with a mass spectrometer shows the element 
to be monoisotopic to 1 part in 50,000. The relative abundances of the molecules 
resulting from electron bombardment of natural phosphorus are given, and the 
existence of the molecule P8 in nature is reported. 


GENERAL 

The most recent tables of isotopes, for example (6), give as a reference for 
the element phosphorus a paper by Aston (1) which dates back to 1920. He 
then reported the element as being apparently monoisotopic. An examination 
of his results leads to the conclusion that at the time he would not have 
observed any isotopes of mass greater than 31 whose abundance was less than 
about 2%, because of the presence of hydrides. His upper limit for isotopes 
of mass less than 31 would have been somewhat better. Since then, no one 
seems to have investigated further, although it was suggested by Naudé (3) 
that a heavier isotope might be present. There was an inviting vacancy at 
mass 33 until the recent discovery by Yaffe and Brown (7) of the radioactive 
P33. The atomic mass of phosphorus as determined by mass spectrographic 
measurements agrees so well with that resulting from other studies that it was 
generally considered improbable that other stable isotopes existed. 


ISOTOPIC ABUNDANCES 
We have recently checked the isotopic constitution of phosphorus, using 
the Laval inflection spectrometer (2). Phosphorus atoms are obtained by 
packing red phosphorus into a filament-furnace which then supplies atoms 
to the ionizing chamber of the spectrometer source. The element is found to 
be indeed monoisotopic to 1 part in 50,000 at least, isotopes 28, 29, 30, 32, 33, 
and 34 being absent to this extent. 
MOLECULAR ABUNDANCES 
Red phosphorus is generally conceded to exist in nature as the molecule P4. 

When evaporated from the hot filament and bombarded by the standard 
70-volt electrons in the spectrometer source, the molecule decomposes. Over 
a moderate range of filament temperatures the abundances of molecular 
masses are in the following proportion: 

Pl (mass 31): 16 

P2 (mass 62): 50 

P3 (mass 93): 6 

P4 (mass 124): 100 

P5 (mass 155): 0O 

P6 (mass 186): 0O 

P7 (mass 217): 0 

P8 (mass 248): 0.5 


1Manuscript received August 19, 1954. 
Contribution from Physics Department, Faculty of Science, Laval University, Quebec, Que. 
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These data indicate that when the molecule splits, it usually splits in half 
under such conditions, the loss of a single atom being a less likely phenomenon, 
as is well known (4, for example). No P5 or P7 molecules were observed, and 
while there was possibly a trace of P6, its abundance would be less than our 
experimental error. The molecule found at mass 248 does not appear to be 
explainable by anything other than P8, although this molecule has not hitherto 
been observed. This experimental evidence for its existence confirms a sug- 
gestion made two years ago by Pauling and Simonetta (5). In discussing the 
question of phosphorus structure, they have proposed that red phosphorus is 
formed from white phosphorus. One of the bonds of the tetrahedron structure 
of the white phosphorus would break, and the open tetrahedrons would then 
join together in long chains, forming red phosphorus. These chains would 
preferably break up into molecules whose mass was a multiple of four, and 
thus P8, P12, etc., might be expected to be found in low abundances. The 
discovery of the P8 fragment confirms this suggestion which was made to 
explain certain anomalies in bond energy measurements. We have not looked 
for P12 (mass 372) but since it would probably be much less abundant than 
P8, it is doubtful whether we would have observed any. 


The author is pleased to acknowledge the financial assistance of the Office 
de Recherches scientifiques of the Province of Quebec in carrying out this work. 
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A POTENTIAL FUNCTION FOR LIQUID HELIUM AT 0°K.! 


By C. F. A. BEAUMONT 


ABSTRACT 


A new potential function for liquid helium is obtained by modifying the 
Margenau potential function and summing over a suggested structure for the 
liquid. The new potential function leads to fair agreement with the first peak 
of the radial distribution curve for liquid helium, with the isothermal compressi- 
bility, and with second virial coefficient data at high temperatures. 


INTRODUCTION 


It is the purpose of this note to determine a potential function between 
two helium atoms in liquid helium at 0° K. Corner and Lennard-Jones (1) 
have given a method for determining a potential function from an experi- 
mental radial distribution curve ; however, the accuracy of the known radial 
distribution curve for liquid helium (4) does not justify the necessary calcu- 
lations. Hildebrand and Wood (3), and Kerr and Lund (6), have given methods 
for determining a potential function of the Lennard-Jones form from radial 
distribution curves and relevant data. This form of the potential function is 
in terms of the maximum interaction energy and collision diameter of the 
atoms; however, it is desirable to obtain a potential function in terms of the 
fundamental forces acting between the atoms. Margenau (9), taking into 
consideration exchange forces and van der Waal’s attractive forces, has de- 
termined the following potential function between two helium atoms: 


[1] o(r) = [770e-*" — 560e-*38" — (1.39/r*) — (3.0/r8)]10-” ergs. 


For purposes of calculation Margenau states that one may use: 
[2] o(r) = [770e-*°" — (1.89/r°) — (3.7/r*)]10—-” ergs. 


If [2] is summed over the structure suggested by Hutchison and Reekie (4) 
one obtains the total potential energy per mole of liquid helium, and a graph 
of this is plotted against the molecular volume in Fig. 1, curve A. The minimum 
value of —96.5 cal./mole occurs at a molecular volume of 16.2 cm.’ as com- 
pared with the accepted value of —14.2 cal./mole (5, pp. 230) at a molecular 
volume of 27.6 cm. (5, p. 207). This indicates that, for liquid helium, some 
additional repulsive potential has not been taken into consideration. 

Yntema and Schneider (12) have developed an empirical potential function 
of the form 

o(r) = aew"”” — (c,/r®) — (¢2/r°) 


to give a good fit with experimental virial data over the range 0°C. to 
1200° C. c; was taken to be 1.24 X 10-” erg-(A)®, as derived by London (7). 
C2, calculated by Margenau (10), corresponding to the above value for ¢,, has 
the value 1.89 X 107 erg-(A)8. a and b were determined by a trial and error 


1Manuscript received September 2, 1954. 
Contribution from the Royal Military College, Kingston, Ontario. 
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Curves A, B, C, and D give the potential energy per mole as a function of the 
molecular volume of liquid helium at 0° K. These curves are obtained by a summation of 
equations [2], [3], [10], and [11], respectively, over the suggested structure (4) for liquid helium. 


=| 
EQUATION [4] A= 0.6406 X10 2 gs, 


e-| 
B=1,.965A 


/ 


-12 
EQUATION C53" =6,C =1,062XI0 ergs 


EQUATION C3] YNTEMA 
AND SCHNEIDER 






EQUATION C23 MARGENAU 


3 4 3 
r in angstroms 


These curves show the potential function between two helium atoms as given by 
Margenau’s theoretical potential function [2], Yntema and Schneider’s empirical potential 
function [3], and the two potential functions [4] and [5], obtained by combining Margenau’s 
potential function with experimental data at 0° K. 
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process to obtain a satisfactory fit with the virial coefficient data. The final 

expression for the potential function was 

[3] o(r) = [1200e~-"/°2!2 — (1.24/r®) — (1.89/r*)]10-" ergs. 

Summing [3] over the suggested liquid helium structure gives a minimum 

value of —55.1 cal./mole at a molecular volume of 19.7 cm.*, curve B, Fig. 1. 
Equations [2] and [3] are shown in Fig. 2 with minima occurring at 2.8 A 

and 3 A respectively. According to the first peak of the radial distribution 

curve for liquid helium (4) the minimum should occur at about 3.44 A. 


DEVELOPMENT 


Equation [2] will be modified by adding to it an additional repulsive term 
either of the form Ae~®’ or of the form C/r”. Equation [2] then becomes 


[4] o(r) = Ae-®" + [770e-*6"" — (1.39/r°) — (3.7/r8)]10-” ergs 
or 
[5] o(r) = C/r™ + [770e-*8" — (1.39/r®) — (3.7/r8)]10-” ergs. 


Considering equation [4] first, it is summed over the suggested structure 
and one obtains 
[6] Wr) = {4.4 (e—% 8909BR + 1.5e7!260BR + Qe} 543BR + Se—!782BR) 
oe [3080(e-* 998% + 1.5e75 7962 oe Qe—7098R + Se—8196F) 
— (15.254/R*) — (42.476/R%)]10-"} N/2 ergs/mole, 


where W(R) is the total potential energy per mole, R* is the volume per atom 
in cubic angstroms, and N is Avogadro’s number. The summation for the 
attractive terms is taken over the 20 nearest neighbors, while for the repulsive 
terms it is sufficient to sum over the first four nearest neighbors. From the 
facts that W(R) = —14.2 cal./mole at R = 3.578 A (molecular volume 
= 27.6 cm.) and dW(R)/dR = 0 at R = 3.578 A one obtains two equations 
for A and B which, on solving, yield the values 
A = 0.6406 X 10-” ergs, B = 1.965 A-". 


Curve C, Fig. 1, is a graph of [6] using these values for A and B. 
At this stage, the isothermal compressibility x7 may be computed from 
[7] 1/xr = V(dR/dV)? (d?W(R)/dR?) 


giving xr = 7.2 X 10-* atm.—' as compared with the extrapolated value of 
12.5 X 10-" atm.—! (2). For a general value of m, equation [5] is not simply 
summed over the crystal lattice to obtain an equation similar to [6]. However, 
summing over the lattice for various particular values of m one may determine 
the corresponding values for C from W(R) = —14.2 cal./mole at R = 3.578 A. 
The correct » is the one giving a minimal value for W(R) at R = 3.578 A, 
and this turns out to be = 6.65. Owing to the difficulties of summation for 
n not an integer, » will be given the value 6 or 7. = 7 yields an isothermal 
compressibility of approximately 5 X 10-* atm~'. For n = 6, C has the value 
1.062 X 10-" ergs, giving an isothermal compressibility of 8.4 * 10-* atm7'. 
Curve D, Fig. 1, shows the total potential per mole for the case m = 6. From 
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Fig. 2 it may be seen that equation [4] has a minimum at r = 3.45 A, while 

equation [5] (with » = 6) has a minimum at r = 3.2 A. Equation [4] is in close 

agreement with the first peak of the radial distribution curve for liquid helium. 
As a further check on the values A = 0.6406 X 10-” ergs, B = 1.965 A-, 

or » = 6, C = 1.062 X 10-” ergs, the second virial coefficients for high 

temperatures may be computed from the classical expression 

[8] B(T) = 2aNJ, (1 — e-# #7) 92 dr, 


where B(T) is the second virial coefficient, N is Avogadro’s number, and & is 
Boltzmann’s constant. A comparison is made with experimental values (11) 
in Table I. 

TABLE I 


COMPARISON OF CALCULATED VALUES OF SECOND VIRIAL 
COEFFICIENTS FOR LIQUID HELIUM 





B(T) in cm.?/mole 





Temp., Measured Equation Equation 
deg. C. values [4] (5] 


0 11.77 12.96 13.06 
300 10.76 11.47 11.05 
600 9.82 10.30 9.97 





One may be tempted to call the difference between the potential energy 
due to summing equation [2] over the lattice and the potential energy resulting 


from summing [4] or [5] over the lattice the zero point energy of liquid helium. 
However, one must note that additivity is assumed throughout the above 
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work. Margenau (9) has shown for helium atoms that additivity is almost 
legitimate for first-order exchange forces and van der Waal’s attractive forces. 
Additivity does not necessarily hold when one considers the zero point energy 
as a repulsive potential energy. Hence curves E and F in Fig. 3, obtained by 
taking the difference between curves A and C and curves A and D respectively 
in Fig. 1, may at best be considered only as a rough approximation to the 
zero point energy. It is interesting to compare these two curves with curve G 
as computed from London’s formula (8) 

[9] E, = h?d/2xm(R — 2-'d)? (R + kd), 


where k = 0.713, R is the mean distance between atoms, and d is the diameter 
of the atoms. In calculating curve G, d is taken as 2.55 A. At a molecular 
volume of 27.6 cm.’ curves E and F give an estimated value of 39.5 cal./mole 
for the zero point energy, while London’s formula gives a value of 41 cal./mole. 


DISCUSSION 
Two possible potential functions for liquid helium have been determined : 
[10] o(r) = [0.6406e—'9*" +.770e—*°9" — (1.39/r°) — (3.7/r*)]10-" ergs 


and 
[11] o(r) = [(1.062/r5) + 770e-*8" — (1.39/r°) — (3.7/r8)]10-" ergs. 


These equations have been constructed to give the correct latent heat and 
molecular volume at 0° K. Owing to the fact that m is integral, equation [11] 
is not quite correct, giving a latent heat of —14.5 cal./mole at a molecular 
volume of 25.5 cm.* With regard to the isothermal compressibility, a sensitive 
quantity, both [10] and [11] yield values which are of the correct order of 
magnitude. Also, from Table I, the agreement with second virial coefficient 
data from 0° C. to 600° C. is within 11% error and the error decreases with 
increasing temperature. Finally the position of the minimum of equations [10] 
and [11] should be about 3.44 A as determined by the first peak of the radial 
distribution curve for liquid helium. The minimum of equation [10] occurs at 
3.45 A and that of equation [11] at 3.2 A. Evidence points slightly in favor of 
equation [10] as being the more acceptable potential function between two 
helium atoms in liquid helium. 


The author wishes to thank Dr. T. S. Hutchison and Dr. D. C. Baird for 
helpful discussions. 
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IONIZATION PROBABILITY CURVES USING AN ELECTRON 
SELECTOR 


RESULTS ON N:*, N*, Xet*} 


By E. M. CLARKE? 


ABSTRACT 


An apparatus consisting of an electron energy selector built into a mass 
spectrometer is described. With it, initial ionization yield curves may be obtained. 
The interpretation of these curves is discussed, and the following new measure- 
ments reported: 

Nitrogen: *II,, state of N2* at 1.0+0.2 ev. above initial ionization potential; 

22° state of N* 2.36 ev. above the 4S® state; 
2P° state of N* 1.4 ev. above the 2D° state; 
energy of dissociation of Ne: 9.80 ev. 


Xenon: _ second ionization potential of Xe: 21.09 ev. 


A. INTRODUCTION 


The measurement of ionization and appearance potentials by electron 
impact methods has been practised for a considerable time (17, 12). Ideally, 
this technique would consist in repeated bombardment of a single atom with a 
single electron, and gradual increase of the energy of the latter until the atom 
was seen to ionize. In practice, one must bombard many atoms with many 
electrons, and so several complicating factors are introduced: 
Experimentally: 

(a) The electrons emitted from a source generally do not have all the same 
energy. For instance, those from a hot filament have a Maxwellian energy 
distribution. 

(6) At the onset of ionization it is necessary to detect very small quantities 
of ions. Thus the threshold of the detector may introduce an experimental 
difficulty similar to that described in (d) below. 

In interpreting the measured curves: 

(c) There may be several ionization processes, to various energy levels. If 
the levels are close together they can produce an effect similar to that due to (a). 

(d) Even for a single process, and monoenergetic electrons, the quantity 
of ionization produced after ionization has set in is not definitely known. 
Wannier (18) has theoretically predicted a near-linear function, which is not 
discontinuous however, but asymptotic to the potential axis, as in Fig. 1. If 
this is so, a precise experimental measurement of the onset of ionization is 
theoretically impossible, unless this function is known in advance. 

To illustrate the effect of these four complicating factors, consider a typical 
ionization probability curve, such as in Fig. 2. Then the point P where the 
curve leaves the axis is quite indeterminate. It depends on the sensitivity of 
the detector, and the theoretical nature of the curve. The curved region R may 

1Manuscript received September 10, 1954. 
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1ON CURRENT 
§ON CURRENT 


ELECTRON ENERGY ELECTRON ENERGY 


Fic. 1. Theoretical ionization yield _Fic. 2. Typical experimental ionization 
according to Wannier. yield. 


be caused by the true nature of the curve (Fig. 1), the energy spread of the 
electrons, or a succession of ionizations to various energy levels, masked by 
the former two. 

Therefore, in order to determine ionization potentials: 

(A) Experimental techniques must be used which effectively reduce the 
electron energy spread, and improve the limit of sensitivity of the detector, 
until these factors no longer mask the various ionization processes. This 
however can never permit the precise determination of the ionization potential 
if the actual ionization curve is an unknown asymptotic function as described 
above. 

(B) So an arbitrary interpretation of the experimental curve such as will 
yield useful information must be made. The worth of any such arbitrary 
method is measured by the success with which it gives results which agree with 
spectroscopic values. 

As far as experimental techniques are concerned, Fox et al. (4) and Morrison 
(14) have recently developed methods which have had considerable success in 
overcoming the electron energy spread, and in enabling them to separate out 
various ionization processes. 

As regards the arbitrary interpretation of the measured curves, this is 
‘currently being done by one of several methods: 

(i) The approximate point at which the curve leaves the axis is chosen 
(vanishing current method). 

(ii) The approximately linear section of the curve is extrapolated to meet 

the axis (linear extrapolation method). 
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(iii) The difference in potential of the ionization curves of two elements is 
plotted as a function of intensity and extrapolated to zero intensity (Warren 
method (19)). 

It is the purpose of this paper to describe a new apparatus for overcoming 
the experimental difficulties mentioned, to discuss the various arbitrary methods 
of interpreting the curves obtained, and to report some new measurements. 

B. APPARATUS 

The apparatus used in this work consisted of a 127° electrostatic velocity 
selector which provided a beam of nearly monoenergetic electrons (9), followed 
by a grid to accelerate these electrons to energies sufficient to ionize the sample. 
The ions formed by this source were separated by an inflection mass spectrom- 
eter of very high luminosity (11) and detected by a 12-stage Allen-type 
Be-Cu multiplier (1, 15), followed by a vibrating reed chart recorder (16). 


Electron Selector 

The electron selector, Fig. 3, was constructed of stainless steel strips set into 
grooves cut in plates of Mycalex K. Electrons from a wolfram filament A 
were accelerated by 3 v. between the filament and the first 0.7 mm. by 4 mm. 
slit of the selector S. The inside curved deflection plate B, of 1 cm. radius, was 
maintained at +4.5 v., while the outside curved deflection plate C, of 1.5 cm. 


9 { 2 CM 


Lp LL 


Ze 


Ms 





a 





Fic. 3. Electron energy selector. 
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radius, could be varied from 0 to +3 v., the usual value for optimum work 
being about +1.3 v. as will be seen. Both slits were adjustable, and at the same 
potential, being mounted on the same strip of stainless steel S. The additional 
energy necessary to cause ionization was given to the electrons by the potential 
difference between the 0.5 mm. by 3 mm. exit slit of the selector and the 
entrance grid of the ionization chamber D. A trap T was provided to measure 
the electron current, and when not in use was grounded to the case D. The trap 
current was found to be independent of the electron energy over a range of 
from 7 to 60 v. 

The Mycalex plates were coated with graphite on the area between the 
deflecting plates to prevent surface charges and to provide a more uniform 
field. A cylindrical magnetic shield 0.5 in. thick and 3 in. long in the form of a 
laminated transformer core was placed around the outside of the spectrometer 
head. 

Gas was admitted to the ionization chamber through a 25 uw quartz leak E£. 
The ions formed were drawn out from D by the fringing field from the grid F, 
and then accelerated by 1000 v. into the spectrometer entrance slits G. Plates 
H and H’ were used for focusing and alignment. The grids on the ionization 
chamber D and the draw-out grid F were of 95% transparency wolfram mesh. 

The ionization voltage was stabilized by cascaded voltage regulator tubes 
VR, and adjusted by a step potentiometer, followed by a potentiometer P 
driven by a motor M, so arranged as to give sweeps of 2, 4, or 8 v. in five 
minutes. The voltages were measured with a volt box and Leeds and Northrup 
potentiometer, calibration marks being keyed onto the motor-driven chart at 
0.20 or 0.50 v. intervals. 

The approximate electron energy distributions furnished by the selector 
were examined by placing a negative bias of 13.5 v. on the trap and measuring 


& TRAP CURRENT x107!° amp 





12,0 13.0 14.0 
ACCELERATION POTENTIAL (VOLTS) 
(Q) 1.1 VOLT ~—O— 
OUTER DEFLECTION PLATE (b)1.2 " —#— 
POTENTIAL (C)1.3 7" —O— 
@)1.4¢ °° —j— 


Fic. 4. Energy distribution furnished by electron selector. 
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the trap current as a function of the acceleration potential. The total current 
curve so obtained was differentiated to obtain the electron energy distribution 
as shown in Fig. 4. 

Here we see the energy distributions furnished by the selector for different 
settings of its outer deflection plate voltage. All of the curves show a sharp 
cutoff on the high-energy side of the distribution, and it is this characteristic 
which is most useful, since it is the long high-energy tail of the Maxwellian 
distribution which causes most of the experimental troubles. The cutoff on 
the low-energy end is due to the relatively sharp cutoff inherent in the Max- 
wellian distribution itself at low energies. Adjusting the outer electrode voltage 
selected different parts of the Maxwellian distribution, and thus more or less 
energy spread was selected. In the figure, total distributions of from 0.8 to 
2.5 ev. are obtained. Thus the combination of selector characteristics and 
Maxwellian distribution of the initial energies results in a smaller energy 
spread than would be expected. 

In addition to this, since the probability of ionization increases with excess 
voltage, the high-energy electrons which appear at the left of the graph are 
more efficient in producing ions, and thus distributions such as Fig. 4(c) and 
(d) are effectively narrower than they appear. It has been found useful in this 
work to take as a measure of the energy distribution the difference in ioniza- 
tion potential as found by linear extrapolation and by the vanishing current 
method. On this basis, distributions as low as 0.2 ev. have been obtained. 
Generally the outer deflection plate was operated at 1.30 v., the trap current 
was about 10-* amp., and the energy distribution about 0.3 ev. 


Mass Spectrometer 

To operate a mass spectrometer with a trap current of 10-§ amp. requires 
special care. The fringing field from an external electrode was used to draw 
out the ions instead of the usual repeller plate, this being found more efficient 
when it was desired to have minimum effect on the electron energy. The grid 
on the draw-out electrode served merely to shield the source from the main 
accelerating potential. The geometric and electrical settings of the focusing 
slits projected ions from a wide angle onto the spectrometer entrance slit. The 
inflection spectrometer was capable of focusing at 12° ion beam onto the detec- 
tor. For stability, the exit slit was opened so that the resolving power was 
about 80 instead of the instrument’s normal 300. 


Detector 

In spite of the precautions mentioned, it was necessary to measure ion 
currents of from 107'§ to 10-'* amp. The Be-Cu multiplier used was sensitized 
in an argon atmosphere, and provided a gain of about 10° when operated at 
about 250 v./plate. Its output of from 10-'* to 10-!° amp. was measured by a 
vibrating reed electrometer with an input resistor of 10° ohms. This combina- 
tion was found capable of detecting single ions, and in the range of ion currents 
actually measured, a 500 picofarad condenser across the electrometer input 
served to smooth out statistical variations in the current. Optimum operation 
of the multiplier was achieved by shielding it from the spectrometer magnet 
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with a 0.5 in. soft iron shield, and positioning a small permanent magnet 
outside the multiplier. These latter precautions improved performance about 


60 times. 
C. INTERPRETATION 


Fig. 5 showing measurements made on xenon and Fig. 6 giving some results 
with NO will serve for the discussion of the interpretation of the experimental 
curves. 

Curves a, 6, and c of Fig. 5 show ionization probability curves of xenon made 
at sensitivities of 1, 1/10, and 1/100. They illustrate the weakness of using the 
vanishing current method of interpreting the curves. Not only are the initial 
points absolutely indeterminate, but they give relative values of ionization 
potential which vary approximately 0.16 ev. for each factor of 10 in the sen- 
sitivity. 

The linear extrapolation method applied to curves 6 and ¢ gives values 
within about 0.02 v. of each other, and is therefore more reliable. However the 
method cannot be used at high sensitivity, such as for curve a, because here 
the line is curved throughout its length, and there is no straight portion to 
extrapolate. In these curves, the electron energy distribution is seen to be 
about 0.3 ev. For this method then, too great sensitivity must not be used if 
good results are to be obtained. Ten runs made on krypton gave in this way a 
precision of measurement of about 0.013 ev. 

Using a Warren plot, the krypton curves gave a precision of measurement of 
about 0.010 ev. However an extra manipulation was necessary to replot the 
curves so as to have the same final slope. 
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Fic. 5. Initial ionization yield of Xe*. Fic. 6. Initial ionization yield of NO*. 
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To summarize, the vanishing current method is poor, the linear extrapola- 
tion method is good if the sensitivity is not too high, and the Warren method 
gives a little better result for some extra work. 

A method not reported to date is the application of Wannier’s theory to the 
curves. As a matter of fact, the linear portions of the curves were rarely 
perfectly straight, but became more so when replotted to give the 1.127 root 
on the basis of Wannier’s prediction that the initial probability should follow 
a 1.127 power law. Using this method, the above runs on krypton gave a 
precision of measurement of 0.009 ev. However this is not conclusive proof of 
Wannier’s thesis, as both the electron distribution and the use of the condenser 
filter might tend to cause a similar curvature. 

So far we have dealt with single processes, where the electron energy distri- 
bution, if small, is of no consequence. Curve 5c shows that the initial approxi- 
mately linear section for xenon is followed by a second curved and linear 
section, indicating that another process of ionization has set in, leaving the 
Xe?* in the ?P 1/2 instead of the *P 3,2 state. Curve 5d corresponding to this second 
process was obtained by subtraction of the linear extrapolation of the first 
straight section from the original combined curve. In five runs, using linear 
extrapolation, the mean difference for the two ground states of Xe+ was 
calculated to be 1.20-++0.05 ev., instead of the spectroscopic value of 1.31 ev. 
This, along with the unduly long tail of the second process, indicates that the 
graphs were slightly curved. However the assumption that they were of the 
1.127 power gave values that were almost as much too high. This figure com- 
pares with the xenon curves obtained by Fox et al. (4) using a pulsed grid 
method to obtain nearly monoenergetic electrons. 

In this case the two processes were clearly separated, but if the separation is 
less than the electron energy spread, or than the initial Wannier curvature, 
then the processes become masked. This is illustrated by argon which also 
has two processes, but separated by only 0.18 ev., so that they merge to give a 
curve similar to 50, instead of 5c. Here linear extrapolation will give an inter- 
mediate value for the ionization potential (12). This is verified in the present 
work: three runs on argon give by linear extrapolation a difference of 3.70 
+.02 ev. compared with xenon, instead of the spectroscopic value of 3.63 ev., 
an error of approximately half the 0.18 ev. in the proper direction. On the other 
hand, when an argon curve assumed to be made up of two processes and a 
xenon curve arbitrarily chosen to have a slope equal to one-half that of the 
argon curve and hence approximately equal to that of the initial process only 
were compared using a Warren plot, the correct difference was obtained within 
0.02 ev. 

The curve for NO, shown in Fig. 6, may be treated similarly. Here, at low 
sensitivity, the line is still curved throughout, and actually fits a 1.7 power law. 
However in this case the curvature is due to the multiplicity of vibrational 
states separated by approximately 0.3 ev., and therefore masked by the 
electron energy spread used. However, as in the case of argon, useful informa- 
tion may be obtained. It is known that the initial 0.5 ev. of the argon curve is 
made up of two processes; and also it is known (2) that there are two processes 
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for NO in the first 0.5 ev. above threshold. Therefore it is interesting to com- 
pare these initial sections of the two curves. When this was done by a Warren 
plot, the ionization potential of NO was found to be 9.24+.04 ev. which corre- 
sponds to the 9.23+.02 value found by Watanabe et al. (20) using a photo- 
ionization technique. Hagstrum’s value of 9.4 ev. (5) approaches the above 
value if his neon curve is reduced to correspond to the slope of only the initial 
processes in NO. 
D. NEW MEASUREMENTS 

Nitrogen: N+ from N» 

Fig. 7 shows a typical curve of N* obtained from Ne. There are two distinct 
breaks in the curve. The first occurs at 2.836+0.07 ev. above the initial ioniza- 
tion potential (mean from 10 curves by the linear extrapolation method) and 







1ON CURRENTS (ARBITRARY UNITS) 


23.00 25.00 ° 
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Fic. 7. Initial ionization yield of N+ from No». 


the second occurs 1.4+0.1 ev. higher. This corresponds very closely to the 
spectroscopic value of 2.38 ev. for the difference in energy between the 4S° 
and 2D° levels of N, and approximately to the additional 1.19 ev. required for 
the ?P° level. From this we may calculate the energy of dissociation of N2 as 
follows. In the process 


N2(X? 23) > Nt+(*P) + N(4S°) 
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(notation as in Hagstrum (5)) we have the energy relation 
Ao = D(N2) + I(N) + E,(N*) ae E.(N) 


where Ay = appearance potential, D = dissociation energy, 


I = ionization potential, E, = electronic excitation energy; 


whence 
D(N2) = Ap — I(N) — E,(Nt) — E,(N). 


The mean value for A» taken from several curves as in Fig. 7 and compared 
with xenon as standard is 24.35+0.1. The ionization potential of N is known 
to be 14.54 (7). E.(N+) and E,(N) are zero for the onset of the process, and 
so we have: 

D(N2) = 24.35 — 14.54 —0 —0 

= 9.81 ev. 

Similarly in the second process 

N2(X? 2}) > N+(@P) + N*(?2D°), 
we have Ao equal to 26.71+0.1 from measured curves as in Fig. 7, and so: 

D(N2) = 26.71 — 14.54 — 0 — 2.38 
= 9.79 ev. 


This result favors the 9.764 ev. spectroscopic value in agreement with 
Douglas (3), and the value 7.4 ev. as reported by Hagstrum (5) would appear 
to be erroneous. Examination of Hagstrum’s paper indicates that he placed 
too much reliance on a value for the heat of formation of NO taken from 
Hanson (6). 

Nitrogen: Nt 

In Fig. 8 we see the initial part of the ionization curve of N.*+ from Ne. The 
curvature from A to C is attributed to the excitation of the vibrational levels 
of the ionized molecule, which are known from spectroscopy (3). Thus if this 
were the only process involved, the curve would continue in a fashion similar 
to that of NO (Fig. 6) as shown by the dotted line. Instead, the curve trans- 
forms into a linear function of the excess voltage along EF. If the curvature 
due to the vibrational levels is ignored for the moment, the portions AB and 
EF of the curve may be extrapolated to meet at D, where another process 
evidently sets in. Because of the nature of the curve, D may not be fixed more 
accurately than 1.0+0.2 ev. above the initial ionization potential. However 
this agrees with the break found by Fox’s method by Frost and McDowell 
(to be published) at 1.20 ev. and attributed by them to the *II, state recently 
identified by Meinel (13) and Dalby and Douglas (2) and situated at 1.115 ev. 
above the 15.576 ev. first ionization level. This interpretation indicates that 
once the electrons have sufficient energy to excite the *II, state, then vibra- 
tional levels are no longer excited. This is not in disagreement with optical 
spectrometry, where the bands attributed to the *II, state have recently been 
found to be quite weak (2). 
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That part of the curve from C to E is attributed to excitation of vibrational 
levels by the electrons in the tail of the distribution of Fig. 4(¢) which have 
not enough energy to excite the *II, state. This effect shifts the apparent 
break from D to E. The curves of Frost and McDowell show a relatively low 
ionization for the first process, and relatively strong ionization for the second 
process. The reverse is true in the present work, and is believed to be due to 
the relatively wide band of electrons exciting a greater number of resonant 
vibrational states at one time. 
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Fic. 8. Initial ionization yield of Ne*. Fic. 9. Initial ionization yield of Xet*. 


Xenon: Xet+ 

Fig. 9 shows a typical ionization curve for Xet*. The line is curved through- 
out, and so the usual methods of interpreting it may not be used. However 
this curvature is not due to a multiplicity of states, as it is known from spec- 
troscopy that the first two excited states of this ion are the *P! and *P°, at 
1.214 and 1.008 v. above the *P? ground state respectively (10). 

Therefore the initial part of the curve should be due to the process involving 
the ground state only. It apparently does not obey a linear law, nor a Wannier 
. function. Arbitrarily, a square root plot, as shown in the figure, gives a straight 
line with a slight break at 1.0+-0.2 v. above threshold, which point seems to 
correspond with the onset of the *P° state. Therefore, taking this initial part 
of the square root plot as representing the ground state process, and com- 
paring it with Xet as a standard, we obtain for the second ionization potential 
of xenon, the value 21.09+0.05 ev. The spectroscopic value is 21.204 ev. (7) 
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and so we may consider the interpretation as possible. Similar approximately 
square-law curves were also found for argon. 

These results are in contradiction with the findings of Hickam e al. (8) 
who have found a correlation with a linear threshold law using the Fox method. 
Following their interpretation, it is also possible that the curve of Fig. 9 is due 
to a linear probability function, plus a large contribution from resonant 
autoionization mechanisms excited by the larger electron energy band used. 

The results on Zet++ and N+ indicate that it is unwise to make a priori 
assumptions of the nature of the probability function; and also that it is 
necessary to use more than one width of electron energy spread to obtain 
maximum information. 
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CALCULATED DIFFRACTION PATTERNS OF 
DIELECTRIC RODS AT CENTIMETRIC WAVELENGTHS! 


By CHARLOTTE FROESE AND JAMES R. WAIT 


ABSTRACT 
Diffraction patterns are calculated for a cylindrical electromagnetic wave 
incident normally on a dielectric cylinder of infinite length whose diameter is 
comparable to the wavelength. Reasonably good agreement is obtained when 
the results are compared with published experimental results. 
INTRODUCTION 

Several experimental investigations of diffraction of electromagnetic waves 
by cylindrical rods whose diameter is comparable to a wavelength have been 
reported recently (1, 2, 7). It was observed that the effects of diffraction around 
the cylinder are very pronounced. In particular, the fields just behind a 
dielectric cylinder exhibit extreme maxima and minima which are difficult to 
explain by any simple physical argument. It is possible that some of the 
effects observed by the above-mentioned investigators might be explained very 
crudely by geometrical optics. In general, however, since the diameter of the 
cylinder is of the order of a wavelength, a complete quantitative description 
can only be found from a solution of the boundary value problem. 

The above-mentioned investigators did carry out calculations based on the 
solution of the boundary value problem for the scalar wave diffraction of a 
perfectly conducting cylinder of infinite length. They found very good agree- 
ment with observed diffraction patterns of brass rods. To the authors’ know- 
ledge no calculations for dielectric cylindrical rods have been published. 
It is the purpose of this note to present such computations and to select the 
parameters to correspond to the experimental conditions. The reader is referred 
to the original papers for details of the observational technique. 


OUTLINE OF THEORY 


The geometry of the diffraction problem is shown in Fig. 1. The source of 
the radiation is taken to be S, the geometrical apex of the horn. SX is the 





Fic. 1. Schematic representation of the dielectric cylinder of radius a with respect to the 
(x, y, 2) and (p, ¢, 2) coordinate systems. 
1Manuscript received August 27, 1954. 
Contribution from the Radio Physics Laboratory of the Defence Research Board, Ottawa, 
Ontario. 
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geometrical axis of symmetry of the horn. The diffracting cylinder of radius a 
is located at O and is coaxial with the z axis of a cartesian coordinate system. 
There are two important cases to investigate, namely when the polarization 
of the incident electric vector is parallel and when it is perpendicular to the 
axis of the cylinder. The parallel case is considered first. 

The distance xo is much greater than x and y in the experiments so that the 
incident wave at O is almost plane, with a slight curvature in the wave front. 
As shown by Kodis (2) a good approximation for the plane z = 0 is to assume 
that the incident wave is cylindrical, which is equivalent to representing the 
excitation by a line source of electric current at S. The exact solution for 
this problem has been given previously by Wait (5) for a homogeneous cylinder 
of arbitrary dielectric constant e, and permeability 4; embedded in a homo- 
geneous medium (¢€o, uo). Employing Watson's (6) notation for Bessel functions 
the field E, at P expressed as a ratio of the incident field Eo at P is* 


(1] Bt = 64 Dy Cyi"H® (kp) Amcos md 
40 m=0 
where 
bo = Bete) by — (x-+x0)), 
k = 2n/X, \ = free space wavelength, 
Co = ie Ci = 2(m ye 0), 
p = (x*+y*)?, — @ = tan“"(y/x), 
and 


bind __ Im (ka) Im(Nka) — Im’ (ka) Im(Nka) 
mn nH (ka) Im’ (Nka) —Ho” (ka) Im(Nka) 


where m = (uoe1/mr€0)?, N= (€1/€0)* (u1/po)?. 





In writing equation [1] in this form, it has been assumed that 
[ a ea STE (1) 

Hy kV (x+x0) +71 io Ho [kxo] im 

Heike +xy) oO" Oe 


which are justified since kx» > 1 and x and y K xo. For purposes of computa- 
tion it is convenient to rewrite equation [1] in the following form 


[2] a = e*— > Cie "sin mm He (kp)cos m¢, 
0 m=0 
where 
Te tan an(ka) — tan am (Nk) 
tan t = tan &, (ba) tan B,»,(ka)—tan am(Nka) ’ 
with Hi = Ho 
and tancd. (ae) = =—Te)y¥ Vial), 





tana, (4) = —x2T— (4)/In(%), 
tan §,,(x) = —xY,,' (x)/ Y,, (x). 


*A time factor exp( —iwt) is employed. 
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When N — i}, corresponding to a perfectly conducting cylinder, it can be 
seen that n», — 6,,(ka). 

When the incident electric vector is perpendicular to the axis of the cylinder, 
as in the experimental work of Keys (1), the ratio of the magnetic field H, 
at P to the incident field Ho at P is identical in form to equation [1] with the 
exception that the coefficient m is replaced by m~! where it occurs in A». 
It is convenient to express the field in the following form for purposes of 
computation: 

[3] Hy _ e”—>> Cit te" "Sin nm Hs (kp)cos m4, 
Ho m=0 

where 

N’tan am (ka) — tan am (Nka) 

N’tan Bm (ka) — tan am (Nka) 


with y: = uo. When N— i! corresponding to the perfectly conducting 
cylinder, it follows that 





tan mm = tan 5m (ka) 


1m* — bm’ (Ra) 


where 
tan 8’(x) = —Jn'(x)/ Yn’ (x). 


In the experimental work of Keys (1), the E, component of the field is meas- 
ured, which is related to H, by 
oe 


¥ few Ox * 


This component of the electric field, expressed as a ratio of the incident field, 
is then written in suitable form for computation as follows: 


[4] = = e*_— ai Curie "sin tm Fn, 
where 
a a m sin m¢ sin oH (kp) 
F, = cosm@ cos ¢H,y, (kp) + ge ee ‘ 


The functions an(x), Bm(x), 5m(x), and 65,’(x) are tabulated extensively by 
Lowan et al. (3). These tables, which are not readily available, greatly simplify 
the computations. 


PRESENTATION OF RESULTS 


Computations from equation [1] have been carried out for selected values 
_ of the parameters. In Fig. 2 the calculated amplitude and phase of E,/Eo are 
plotted against ky for a polystyrene cylindrical rod for x» = 600 cm., A = 1.25 
cm., @ = 1.27 cm., and with the relative dielectric constant N? taken to be 
2.56 to correspond with the experimental conditions of Kodis (2). The experi- 
mental curves are shown by the broken curves which correspond quite closely 
with the calculated values except for small values of ky. The discrepancy here, 
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Fic. 2. The amplitude and phase of E./Eo for a polystyrene cylinder whose radius is 
approximately equal to a wavelength. 


which is not more than 13%, can probably be attributed to the ohmic loss 
in the polystyrene which has been neglected in the calculation. Actually, 
the dielectric constant N? according to published values (4) would be 
2.56(1 + 70.0012) at a wavelength of 1.25 cm. It is also possible that the inter- 
action between the probe and the cylinder might account for this disagreement. 

Similar computations from equation [1] have also been carried out for a 
lucite rod (N? = 2.56, a = 1.22 cm.) and a tenite rod (N? = 3.20, a = 1.27 
cm.), both for x-+x»9 = 615 cm. and \ = 3.2 cm. These results are plotted in 
Figs. 3 and 4. Experimental values of McLay are shown in Fig. 3 by the 
broken curves which correspond very closely to the computed curves. It is 
also interesting to compare the computed patterns for the lossless tenite 
cylinder in Fig. 4 with Wiles’ and McLay’s (7) published experimental results 
for a hard rubber cylinder. All conditions for the two cases are nearly identical 
except that the rubber has a loss tangent of about 0.10, which has the effect of 
greatly reducing the maxima of intensity just behind the cylinder. The other 
aspects of the curve, however, are very similar, which would substantiate our 
explanation that the focusing action of the cylinder is seriously affected by the 
resistive loss in the dielectric. 


tProf. McLay has recently repeated the experiment with a tenite cylinder and the results are 
within a few per cent of the calculated values shown in Fig. 4 
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Fic. 3. The amplitude of (E./Eo)? for a lucite cylinder whose diameter is comparable to a 
wavelength. 


The complexity of the calculations for the perpendicular polarization is 
somewhat greater. The computations from equation [4] for x+x) = 535 cm., 
A = 3.2 cm., a = 1.22 cm., N? = 2.56, are illustrated in Fig. 5 for several 
values of x. Experimental values of Keys (1) are shown by the broken curves. 
. The agreement is reasonably good except for points in the plane x = 0. 
It is believed that the discrepancy here can be partly attributed to the spurious 
response of the detector probe to the E, component of the field. In the plane 
x = 0, the E, component of the scattered field from the cylinder is much larger 
than the E, component so that slight misalignments or asymmetry of the 
probe will produce a response which is not indicative of Fy. 








CANADIAN JOURNAL OF PHYSICS. VOL. 32 







DIFFRACTION BY TENITE CYLINDER 
(it POLARIZATION) 


. HAA AACA 
At el dd 
a ee 
LAIN NIAN ALT | | 
CNP 
dalle cs 
















xX* -67Cm. 


INA ANAL 
TV YY 
Uo 
ARREARS 
OT eset en | 
IAAT ALLL 
SO YP 
USEC asetent 
See 
CIA CALALAGA 
CTW EY 
Se eae 



























Dee 
y ccm) 






Fic. 4. The amplitude of (E,/Eo)? for a tenite cylinder whose diameter is comparable to a 
wavelength. 











CONCLUSION 


These results show that a predominant characteristic of the diffraction 
pattern of a dielectric rod is the peaks of intensity near the center of the 
geometrical shadow just behind the cylinder. This effect is apparently present 
when the incident wave is either parallel or perpendicular to the axis of the 
cylinder. These peaks seem to be generally slightly lower in the experimental 
patterns, which is probably a result of resistive loss in the dielectric. For the 
case of lucite and polystyrene the dielectric ‘‘loss tangent” is never greater 
than about 0.005 at centimeter wavelengths and the assumption that the 
dielectric constant is real would not be expected to lead to serious error. 
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Fic. 5. The amplitude of (E,/Eo)? for a lucite cylinder with the polarization of the incident 
electric vector transverse to the axis of the cylinder. 


In all the diffraction patterns it is noted that intensity maxima and minima 
vary in a very regular manner in the region away from the geometrical shadow. 
Furthermore, the location of these maxima and minima is only slightly 
dependent on the dielectric constant of the cylinder. 


ACKNOWLEDGMENT 


Appreciation is expressed to Prof. A. B. McLay who communicated to us 
the experimental results shown in Fig. 3. Helpful suggestions have been given 
by J. H. Crysdale and J. E. Keys. 


REFERENCES 


1. Keys, J. E. M.Sc. Thesis, Dept. of Physics, McMaster University, Hamilton, Ont. 1953. 
To be published. 

2. Kopis, R. D. J. Appl. Phys. 23: 249. 1952. 

3. Lowan, A. N., Morse, P. M., Fesuspacu, H., and Lax, M. Appl. Math. Panel Rept. No. 
62.1 R. Mass. Inst. Technol. Feb. 1945. 


. Tables of dielectric materials. Vols. I-III. Laboratory for Insulation Research of the 
Massachusetts Institute of Technology, Cambridge, Mass. June 1948. 


. Watson, G. N. Theory of Bessel functions. 2nd ed. Cambridge University Press, London. 


1944. 


4 
5. Wait, J. R. Geophysics, 17: 378. 1952. 
6 
7 


. WiLEs, S. T. and McLay, A. B. Can. J. Phys. 32: 372. 1954. 





THE EFFECT OF HIGH LOADING ON THE KINETIC FRICTION 
OF ICE! 


By C. D. NIVEN 


ABSTRACT 

The kinetic friction of stainless steel on ice decreased very markedly at high 
loadings of 30 to 50 kgm. per sq. cm. This constituted a violation of Amontons’ 
Law. At low loading the material of a slider is of far greater importance than at 
high loading, which suggests that in the explanation of the friction of ice at least 
two important factors are involved, which are related respectively to (a) pressure 
and (b) adhesion. 

INTRODUCTION 

Although the low friction of ice just under the freezing point must have 
been observed by man from the earliest times, yet the explanation of the 
phenomenon has been only partially given. Physicists and engineers have been 
for years content with the theory that pressure tended to melt the ice and so 
created a film of water to furnish lubrication. More recently the suggestion 
has been put forward that the liquid for lubrication is provided by frictional 
heat (3) but this at once raises the question as to why certain good metallic 
conductors of heat have a low friction on ice. 

The explanation of this extremely commonplace phenomenon is certainly 
not obvious and needs to be completely reconsidered in the light of modern 
developments in the theory of friction (1, 5, 6). The lack of information at 
high loading is surprising. Interest seems to have centered mainly around the 
ski for personnel or aircraft rather than around the sleigh runner. A report, 
however, by Wakefield (8) on logging sleighs sliding on packed snow, published 
in 1938, indicated that the friction at high loading decreased. In 1939 Bowden 
and Hughes (3) using a turntable showed that the friction of ice dropped when 
the load on the slider was increased above ca. 300 gm. Only loads under 600 gm. 
were used and the area of apparent contact was not stated. The data, none- 
theless, supported Wakefield’s work and to some extent fitted in with the field 
observation that after a certain load had been put on a runner drag increased 
little as the load was further increased. 

Altogether there was a distinct suggestion that insofar as the kinetic friction 
of ice was concerned Amontons’ Law was not always true at high loading. The 
work described below deals with the investigation of this point. 

APPARATUS 

A photograph of the turntable which was used by the writer is shown in 
Fig. 1. It consisted of a 36 in. diameter disk of 3/16 in. boiler plate on the edge 
of which was fixed a circular bath for containing the water to be frozen for the 
track. The boiler plate and the track were mounted at the center on a vertical 
spindle protruding from a gear box. A d-c. motor, with series resistance, was 
used to drive the turntable. 


1Manuscript received August 31, 1954. 
Contribution from the Division of Physics, National Research Council, Ottawa, Canada. 
Issued as N.R.C. No. 3487. 
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Fic. 1. Photograph of turntable, showing tool-holder for machining surface of ice, lever arm 
for putting load on slider, and gravity balance for measuring the drag. Since the circular trough 
was empty at the time of photographing, the arm is not resting on the slider; the latter can be 
seen in the trough near the arm. 


From a point above the center an arm or lever was pivoted so that it could 
be raised above the track when the ice was being machined, and lowered on 
to the sample or slider during a test. The pivot was somewhat like a universal 
joint because provision was made to allow the arm to swing horizontally 
through a small angle. To reduce friction against motion in this plane ball 
bearings were used. The sample supported the arm at a point 15 in. from the 
center or pivot and a hook was provided 25 in. farther out on the arm for 
carrying weights. A 60 lb. weight could by this means put a load of 193 Ib. on 
the slider; on a sample of area 2 in. X 1 in. this approximated 100 Ib./sq. in. 
To support the table under the sample a wheel had to be mounted under the 
boiler plate. 

The track was 6 in. wide and 3 in. deep and after the ice was frozen it had 
to be machined. To do this a tool holder was mounted on the bench of the 
turntable in such a way that it could be swung into place while the ice was 
being surfaced, and swung out of place when the experiments were in progress. 
Fig. 1 shows the tool holder swung out of place and the arm down. 

The balance was a simple gravity balance (see Fig. 2). The sensitivity far 
exceeded the sensitivity required to measure the “‘stick-slip’’ drag of the 
friction test. The speed of rotation of the turntable was kept down to about 
one revolution in three seconds; as the diameter of the circular track was 
33.5 in. this worked out to about 89 cm. per sec. Since 1 mile per hour is 
approximately 45 cm. per sec. the speed was considerably lower than the 
speed which would normally occur in transportation; therefore the drag values 
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Fic. 2. Gravity balance designed to weigh a horizontal drag. 


can be considered higher than would be encountered in practice on similar 
ice surfaces and in the same atmospheric conditions. On account of the fact 
that the slider was running over the same track again and again, speeds were 
kept low intentionally, in order to prevent too much heat from being put into 
the surface of the ice; obviously the work done by the motor on the slider must 
for the most part go into the ice as heat. A 3-in. steel ball resting in a block of 
brass permitted the slider to adjust to any angle in any direction under the 
lever carrying the weights. 

The cold room in which the tests were made could be controlled to +1°F. 

Twenty pounds were as a rule added at a time to the hook at the end of the 
lever and with each addition the drag was balanced with weights on the 
balance pan. Since the actual load on the slider was the thrust of the arm 
itself plus 22 times the weight on the hook, 0, 20, 40, and 60 Ib. on the hook 
corresponded to loads on the slider of 14.92, 39.13, 63.35, and 87.56 kgm. 
These figures are the abscissae of points on the drag vs. load charts. Values of 
the friction constant have not been worked out because the significance of a 
constant rather disappears if it is not a constant at all. 
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RESULTS 


Results on stainless steel, bakelite, and teflon (polytetrafluoroethylene) are 
shown in Figs. 3, 4, 5. These graphs represent drag plotted against load (i) on 
a sample of apparent area 13 sq. cm.—indicated by full lines—and (ii) on a 
sample of apparent area 1.6 sq. cm.—indicated by broken lines. The larger 
apparent area corresponds to the area of a piece of material 2 in. X 1 in. while 
the smaller apparent area corresponds to the area of two runners on a miniature 
sleigh each 1 in. long and ¢ in. wide. All drag vs. load graphs are plotted to the 
same scale so that direct comparisons can be made. 


STAINLESS STEEL 


APPARENT CONTACT AREA 13 SQ, CM. 
“ “ * 6 SQ. C~.— —— 


ORAG - GRAMS 


— ee 
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Fic. 3. Drag vs. load curves for stainless steel. 


The measurements were made at 29°F. and 1°F. The samples all slid on the 
ice fairly smoothly except the larger sample of stainless steel at 1°F. In this 
instance the stick-slip motion, so well known in the measurement of friction, 
became very troublesome and the drag value was most indefinite. A question 
mark against the point on the diagram draws attention to this. Only the 
lowest loading was measured because the jerky motion was so hard on the 
apparatus. 

The graphs are not plotted from averages. Each curve is plotted from a single 
set of readings which is representative of the behavior of the material at that 
temperature. Friction measurements are notoriously difficult to repeat to a 
high degree of accuracy, and with ice slight changes in temperature can cause 
a change in the absolute friction value. It was found however on examining the 
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Fic. 4. Drag vs. load curves for bakelite. 


TEFLON 
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Fic. 5. Drag vs. load curves for teflon. 
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experimental data that these slight variations were quite unimportant in 
drawing several conclusions. These can be summed up as follows: 

(1) All sets of readings showed very marked curving of the graphs for 
samples of small apparent area of contact at 1°F. but little or no curving for 
the samples of larger apparent area of contact at that temperature. 

(2) All sets of readings for teflon and stainless steel samples of small 
apparent area of contact showed an almost constant drag at 29°F. 

(3) The readings consistently showed that the bakelite samples whether of 
small or of large apparent area of contact had a higher drag value at 29°F. than 
the corresponding teflon samples. At 1°F. the accuracy of the work does not 
however permit discrimination between these two materials. 

(4) At both 29°F. and 1°F. the effect of apparent area of contact on the 
kinetic friction is far more marked for stainless steel than for either bakelite 
or teflon. 

(5) Considering only samples of large apparent area of contact, the increase 
of kinetic friction with lowering of temperature is far greater for stainless steel 
than for either bakelite or teflon; on the other hand, considering only the 
samples of small apparent area of contact, while there is an increase in kinetic 
friction with lowering of temperature all three materials are comparable at 
1°F. with a distinct indication that at very high loadings stainless steel may 
slide the best of the three materials. 


DISCUSSION OF RESULTS 


It is evident that under certain circumstances depending upon temperature 
and loading, the field observation that drag for steel is not increased much 
with increase of load at high loadings is essentially correct. In other words 
Amontons’ Law fails at high loading on ice under certain circumstances. Such a 
conclusion for kinetic friction contrasts with Bowden’s (2) conclusion on 
static friction. Teflon, which has come on the market only recently, does not 
show much improvement over bakelite around 0°F., but it is definitely better 
than bakelite just below freezing; bakelite has been advocated by Klein (7) 
for aeroplane shoeing before teflon was discovered, but teflon seems to have a 
possible advantage worth investigating under more realistic conditions than 
the laboratory permits. 

The curving towards the load axis at high loading, which is most clearly 
shown in the case of stainless steel around 0°F., taken in conjunction with the 
fact that stainless steel contrasts severely with the two plastics in certain 
other aspects of the work, suggests the following micromechanism for the low 
friction of ice. After the molecules which constitute the extreme tip of some 
asperity on the ice surface have been forced to start moving by the pressure 
and sheering forces during the sliding process, they have become ‘“‘agitated”’ 
or in other words “‘heated’’; they have received a part of the energy which the 
motor has wasted in dragging something along the level. Were this energy 
wastage great enough one would expect melting, and using a rubber slider and 
a temperature just below freezing this melting of the ice surface is quite easy to 
observe; on the other hand using a hard slider and a lower temperature, it is 
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not easy to detect actual liquid. However it is not really of importance whether 
liquid can be detected either visually or by means of a chemical indicator. The 
important thing is that the molecules on the tips of the asperities have become 
agitated both on the ice surface and on the slider surface; and when molecules 
are agitated they are conditioned—so to speak—for forming a bond if bonding 
is at all feasible. No one would argue this point if one were discussing heating 
reagents in a test tube over a bunsen burner or welding platinum wire, but it 
is not so obvious when the agitation has arisen mechanically through dragging 
one substance over another. With ordinary solid materials such as constantan 
sliding on mild steel, temperatures as high as 1000°C. have been observed by 
Bowden and Ridler (4) during the friction process. The tendency to form 
liquid when one substance slides over the other must then normally assist the 
formation of an adhesive bond rather than assist in the sliding. 

Ice however is not a normal substance and there is, in consequence, a special 
point to consider. If a peak on the ice surface is approaching an inverted peak 
on the slider surface, when they actually make contact there is a sudden 
increase in pressure. How great this is will depend on the load and on the 
velocity of approach. Ice has the peculiar property of tending to melt rather 
than freeze with increase of pressure, so that agitated molecules which are 
available for bonding are discouraged from getting into the solid state during 
the collision of two peaks although this is just the very time that any normal 
substance would find it most easy to take advantage of the possibility of form- 
ing a bond. 

If the temperature is very cold this tendency to melt rather than freeze 
when two peaks collide is too small to consider in comparison with the tendency 
of the surrounding cold to stop agitation and so form a bond. Thus the friction 
of very cold ice would be expected to differ little from that of any other solid. 
Those materials which are hydrophobic would not bond easily with ice even 
when the molecules were agitated, so that there would be no special reason why 
increase of pressure at peaks at the time of collision should help prevent some- 
thing that can only occur with difficulty anyhow. 

Weyl (9) considers ice to be covered with a fluid surface layer hundreds of 
molecules thick. The idea of this thick layer of distorted water molecules would 
fit in much better with the observed friction facts at high loading if this layer 
were not a fluid but a solid which like ordinary ice contracted on turning to a 
liquid. Any theory which fails to make use of this basic and almost unique 
property of ice in explaining its low friction seems questionable. 


CONCLUSIONS 


1. Amontons’ Law is not as a rule true for substances sliding on ice when 
the loading per unit area of apparent contact is high. 

2. There appear to be two entirely different factors at work in determining 
the friction of a substance on ice, firstly a general friction factor which depends 
on the loading and on the temperature and which applies to all hard materials 
sliding on ice, and secondly a special friction factor which depends on parti- 
cular qualities of the sliding material. How these two factors combine to give 
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the resulting friction is not clear but they are certainly not additive. The work 
suggests that the second factor is relatively less important than the first at 
very high loadings per unit area. 
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MOTION OF CLOUDS OF ABNORMAL IONIZATION IN THE 
AURORAL AND POLAR REGIONS! 


By E. L. Hacc anp G. H. HANSON? 


ABSTRACT 


Examination of records taken in rapid succession at ionospheric stations in 
Northern Canada has revealed three distinctive types of echoes corresponding 
to clouds of ionization in motion. One type is ascribed to clouds of sporadic E 
moving horizontally over the station. The second type is believed to be due to 
clouds of ionization descending vertically from the F region to the E region, 
while the third type appears to be from clouds moving at extremely high velocities 
in the E region. Reflections of the second and third types have been seen only at 
stations in or north of the auroral zone. The characteristics of the echoes, the 
velocities of the clouds, and other associated phenomena are discussed. 


I. INTRODUCTION 


Canadian ionospheric stations are provided with automatic recording 
equipment which can obtain complete h’-f records of the ionosphere as often 
as 10 times per minute. For several years such records have been taken for 
24 hr. every two weeks at Baker Lake (64.3°N., 96.0°W.) and at Fort Chimo 
(58.1°N., 68.3°W.). A smaller number are available from Resolute Bay 
(74.7°N., 94.9°W.). Unfortunately, the interference created by the operation 
of the equipment has prevented the taking of continuous runs at the other 
northern station at Churchill (58.8°N., 94.2°W.) and has curtailed the pro- 
gram at Resolute Bay. The h’-f sweeps are photographed on 16 mm. film in 
such a way that when the film is run through a motion-picture projector, a 
speeded-up picture of the changes in the ionosphere is revealed. Many short- 
lived reflections and rapid changes are observed which would be missed 
completely by the routine sweeps which are taken at 15-min. intervals. 

Most of the film from these continuous runs has been carefully examined 
recently in order to study certain unusual types of echoes which had been 
observed earlier. These echoes have been divided into three distinct types. 
The characteristics of each are given in succeeding sections. 

2. VELOCITIES OF CLOUDS OF SPORADIC E 

When viewing continuous runs, traces identical in appearance with sporadic 
E often appear at virtual heights of 200 to 300 km. and then decrease in 
height in a regular manner until they reach the 100 km. level. The traces 
often remain at this level for some time and then increase in height again 
until they disappear. They will be referred to as type I traces. A typical 
graph of height as a function of time is shown in Fig. 1. It is evident that 
type I traces can be interpreted as being due to reflections from a cloud of 
sporadic E moving horizontally at a height of about 100 km., the minimum 

1Manuscript received June 16, 1954. 

Contribution from the Defence Research Telecommunications Establishment (Radio Physics 
Laboratory), Defence Research Board, Ottawa, Canada. Work carried out under DRB Project 
No. D48-97-52-01. 
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Fic. 1. The change in virtual height of a sporadic E cloud as it approaches the station. 


virtual height occurring when the cloud is either directly over the station or 
is at a point of closest approach. This phenomenon has been reported pre- 
viously by Meek (3), who found 40 examples in film from Baker Lake. Of his 
40 examples, 11 were of clouds which appeared to belong to the E region, 
while the others were believed to be higher. 

A total of 51 more occurrences of type I clouds have been observed in this 
study. Of this number, 26 are from Resolute Bay, 15 are from Fort Chimo, 
and 10 are from Baker Lake. Graphs of height vs. time have been drawn 
for each occurrence and the velocities of each have been computed. In cal- 
culating the velocities it was assumed that the cloud was moving horizontally 
at constant speed at a height equal to the lowest height observed, or, if the 
minimum height observed was above 120 km., that the true height was 120 km., 
but that the cloud passed some distance to the side of the station. The dis- 
tribution of velocities for each station is shown in Fig. 2. The median velocity 
is about 800 km. per hour with a few being three times as great. There are 
no significant differences in the velocities measured at the three stations. 

The number of occurrences as a function of time of day is shown in Fig. 3. 
There is no favored time of day at Resolute Bay, but at Baker Lake and Fort 
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Fic. 2. The distribution of velocities of type I clouds. 
Fic. 3. Diurnal distribution of type I clouds. 
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Chimo most type I clouds were observed in the evening. This is in agreement 
with observations of sporadic E made regularly at these stations. At Baker 
Lake and Fort Chimo, both in the auroral zone, most sporadic E occurs at 
night, whereas at Resolute Bay there is no pronounced diurnal variation. 
There were not enough observations to determine whether or not there was a 
diurnal variation of velocity. However, a greater range of velocities occurred 
between midnight and 8 a.m. and all the very high velocities occurred in the 
early morning. 

The velocity of motion of the sporadic E cloud should be that of the air at the 
same height. It is interesting to compare the above results with measurements 
of ionospheric winds made by other methods. Chapman (1) has measured 
winds in the normal E region at Ottawa and Montreal, using the fading of 
radio echoes at spaced receivers. He found the velocities to range from 70 
km. to 1080 km. per hour, with a median value of approximately 290 km. per 
hour. This is somewhat less than was found in this work where the velocities 
range from 270 to 3500 km. per hour and the median value is approximately 
800 km. per hour. It is not too surprising to find these differences, however, 
for the observations were not made at the same location, nor were the reflec- 
tions necessarily from the same height in the ionosphere. 


3. VERTICAL MOTIONS FROM THE F REGION TO THE E REGION 


Another type of moving cloud, hereafter referred to as type II, is illustrated 
by the series of photographs in Fig. 4. These photographs were taken 15 sec. 
apart and show the sudden appearance of a reflecting ‘patch’ below the F2 
layer. The patch moves rapidly downward to a minimum virtual height of 
180 km., and then rises to disappear in the F2 layer. A height vs. time plot of 
the cloud in Fig. 4 is shown in Fig. 5. 

The characteristics of type II reflections are as follows. The trace first 
appears below the F1 or F2 trace and has the same curvature as the F traces. 
There is evidence of retardation at a critical frequency and of splitting into 
o and x branches. As the virtual height decreases, the highest frequency in- 
creases and the trace becomes flatter; when the height increases, the reverse 
changes take place. There are always some spread echoes and fuzziness in con- 
nection with the reflections, and occasionally the height fluctuates for a few 
moments at or near the minimum height. There is no blanketing of the F 
region and no significant effect on the regular F layers. The trace never appears 
at a frequency lower than the critical frequency of the E layer. 

The above characteristics indicate that the cloud of ions causing the reflec- 
tions is moving vertically from above the F2 layer into the E region. If such 
is the case, the cause of the phenomenon is hard to imagine. It is as if a 
stream of some ionizing radiation were suddenly turned on from outside the 
earth, the radiation becoming more intense and penetrating for some minutes 
and then gradually tapering off. 

The velocities of each patch have been computed on the assumption that the 
motion is vertical. The distribution of velocities for the three stations where 
the phenomenon has been seen is shown in Fig. 6. The average downward 
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Fic. 5. The minimum virtual height of the cloud of Fig. 4 as a function of time. 
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Fic. 6. The distribution of velocities of type II clouds. 


N@ OF OCCURRENCES 


velocity is about 7200 km. per hour although in one case a velocity three times 
as great was observed. There does not appear to be any difference among the 
stations although too few examples are available to be certain. 

One significant point that has been noted is that the minimum height 
reached is lower at Fort Chimo and Baker Lake than at Resolute Bay. At the 
former stations it is usual for the minimum height to be 120 to 130 km. while 
at Resolute Bay it is about 150 km. When the ionization does reach to the E 
region, the critical frequency of the normal E£ layer increases slightly, or if 
the event occurs at night, a night E layer (2) forms. There have not been 
enough examples found as yet to determine whether type II is more common 
in the daytime or at night. 

It must be emphasized that type II traces occur only at Fort Chimo, 
Baker Lake, and Resolute Bay. Many continuous runs are available from 
Prince Rupert and Winnipeg, far to the south, but no occurrence of type II 
has been seen there. With respect to the auroral zone, Prince Rupert, Winni- 
peg, and Resolute Bay are nearly equivalent, Resolute Bay being as far to 
the interior as the other stations are to the exterior. 

It should be noted that the 28 occurrences which have been scaled do 
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not give a true indication of the frequency of occurrence. On many other 
occasions, similar traces were seen but could not be accurately scaled because 
of spread echoes, lack of height markers on the film, or various equipment 
troubles. 


4. EXTREMELY RAPID MOTION IN THE E REGION 


A third type of abnormal ionospheric reflection has been observed at Fort 
Chimo and Baker Lake. Type III echoes are in many ways similar to type I 
echoes, except that the apparent velocity is much higher. As well, absorption 
is invariably much increased after the trace reaches the E region, and black- 
outs occur as an aftermath. The appearance of the reflections is shown in the 
series of photographs in Fig. 7, while the rate of change of virtual height can 
be seen from the graph in Fig. 8. The distribution of the velocities for Fort 
Chimo and Baker Lake is shown in Fig. 9. 
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Fic. 8. The virtual height of the type III reflections of Fig. 5 as a function of time. 
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Fic. 9. The distribution of velocities of type III clouds. 


The type III traces usually appear first at virtual heights of from 300 to 
500 km. at a frequency of from 6 to 10 Mc. or more. As the height decreases 
the frequency range of the reflections decreases, although not regularly, 
until at a height of 100 km. the trace usually extends 2 or 3 Mc. beyond the 
foE, that is to about 6 Mc. In some cases the trace has descended to a height 
of 85 km. 

From the appearance of the reflections and from the fact that no retardation 
has been seen near the F2 layer critical frequencies, it would seem that the 
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ion cloud is moving horizontally in the E region. On this assumption, the 
velocities have been computed. The average velocity is about 4500 km. per 
hour but in two cases at Baker Lake the velocity was nearly 15,000 km. per 
hour or 4 km. per second. There were several cases where it was felt that even 
higher velocities occurred, but the trace appeared on only a few frames and so 
was not scaled. 

On most of the occasions when type III reflections were seen, a great 
increase in absorption, resulting in a total blackout, occurred when the trace 
reached a height of 100 km. or less. When a total blackout did not occur, the 
absorption increased (as seen from the minimum frequency of reflection) 
and the E region appeared very turbulent. Often a night E layer (2) appeared, 
and after the absorption returned to normal, the E region remained turbulent. 
The occurrence of blackouts or night E prevented observations of the trace 
after it reached its minimum height. 

Type III reflections have been seen only at Fort Chimo and Baker Lake, 
in the auroral zone. The suspicion that the echoes are connected with the 
aurora is supported by Fig. 10 where the diurnal distribution is shown. 
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| TYPE II - BAKER LAKE 
0 2 8 So 8 & 6B 6 He S&S & 


LOCAL TIME 
Fic. 10. Diurnal distribution of type III echoes. 
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All occurrences were between 4 p.m. and 5 a.m., roughly the same time as the 
aurora is observed. 

One is tempted to suggest that type III echoes occur when a stream of 
charged particles enters the ionosphere and the point of contact of the stream 
with the E region moves horizontally. However, the velocities are very much 
higher than those of auroral structures which last for periods of several 
minutes, and it is hard to see how the echoes could be explained on this basis. 
It is probable that a satisfactory explanation of type III and type II echoes 
cannot be made unless simultaneous measurements of magnetic and auroral 
characteristics at the same station are available. 


5. CONCLUSIONS 


Type I echoes appear to arise from clouds of sporadic E ionization drifting 
with the winds in the E region of the ionosphere. The wind velocities measured 
on this assumption vary from 235 to 3500 km. per hour with a median value 
of 800 km. per hour. This is in general agreement with observations by other 
authors, having regard to the differing places of observation. 

The echoes designated as type II and type III cannot be explained as 
simply as those of type I. Although type I and type III echoes have some 
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similarities, the extremely high velocities of type III echoes (up to 1500 km. 
per hour) make it most unlikely that the apparent motion is due to an iono- 
spheric wind. Since type III reflections have been observed only at stations 
very close to the auroral zone maximum, it is quite possible that the effect 
is due to sweeping of the auroral ionizing agent. Type II echoes, while not as 
strictly limited to the auroral zone, also appear to relate to an auroral ionizing 
agent, and can be explained as due to a stream of ionizing radiation which 
has suddenly been turned on. The apparent downward velocity measured on 
this premise is about 7200 km. per hour. 
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CALCULATION OF DEPOLARIZATION FACTORS IN THE RAMAN 
SPECTRA OF POLYATOMIC MOLECULES! 


By SALVADOR M. FERIGLE AND ALFONS WEBER 


ABSTRACT 


General expressions for the derived polarizability tensor are deduced for 
non-cyclic molecules having bonds with cylindrical symmetry by assuming the 
additivity of bond polarizabilities as suggested by Wolkenstein. Full use is made 
of the Wilson FG matrix technique. As an example, the expressions for the 
trace and the anistropy for the CX.Y molecular type are obtained explicitly. 


INTRODUCTION 


The calculation of intensities and depolarization factors in the Raman 
spectra of polyatomic molecules has been given an increasing amount of 
attention in the last few years. The theory—as it stands today— assumes the 
additivity of bond polarizabilities, and this assumption, together with a few 
complementary assumptions about the dependence of the bond polarizabilities 
on the internal coordinates, was made the basis of the method developed by 
Wolkenstein (10). The same basic ideas have been used by Long (6), who 
used, however, different mathematical techniques more in line with the FG 
matrix method developed by Wilson (9) for the calculation of vibrational 
frequencies of polyatomic molecules. In this paper we attempt to make a more 
complete use of Wilson’s formulas and to give a more explicit formulation 
for the elements of the derived polarizability tensor than was given in previous 
papers. 

A set of rectangular coordinate axes xyz is fixed in the molecule and the 
additivity of bond sei and their derivatives is expressed as 


[1] me a”) 


and 
(n) 
(2] whe 


where u and v can take the values x, y, 2; refers to the mth bond in the 
molecule; the summation is carried over all the bonds in the molecule; and 
Q is the kth normal coordinate. 

The depolarization factor p, of the kth fundamental vibration is given by 
(8) 


(3] _ __ 6B," 
Pe 5Ay TBS 


where A, and B, are respectively the trace and the anisotropy of the derived 
polarizability tensor da,,/dQ*, that is, 
1Manuscript received July 22, 1954. 
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Dvn 
Ar= Lage 


a a Se ) +S y) + Gat) 
(9) 


For the special case of a symmetric top molecule, a,;; = a, ¥ a,, and, for 
totally symmetric vibrations, 


[5] By = A;,—30a;;/dQ*. 


The principal axes of the bond polarizability are assumed to be both in the 
equilibrium position and during vibration to be directed along and perpendi- 
cular to the bond, and therefore 


[6] a = SY eamcos(nmu)cos(nmv), 














where m takes the value 1 for the component along the bond and the values 
2 and 3 for the axes perpendicular to the bond, and cos(mmu) is the direction 
cosine of the mm-axis with the u-axis. For bonds with cylindrical symmetry 
Qn2 = Q,3, and Eq. [6] gives, from the orthonormal properties of the direction 
cosines, 


a? = anycos(n1u)cos(n1v) +en2[cos(n2u)cos(n2v) +cos(n3u)cos(n3v) | 


= (cn1—,2)c0s(21u)cos(n1v) +an25u0; 


where 6,, is the Kronecker delta symbol. If we write (mu) for (mlz) since it is 
simply the angle between the mth bond and the w-axis, the last equation can 
be written as 


[7] as” = (ani —On2)COS(nU)cos (nv) +en28ue- 


Wolkenstein assumed that a», depends only on the stretching of the mth 
bond, that is, it is unaffected by the stretching of other bonds or changes in 
interbond angles. Then, the molecular polarizabilities depend on the bond 
stretchings through a,» and on the interbond bendings through the direction 
cosines. As clearly stated by Long, Wolkenstein’s assumptions are not neces- 
sary and a similar treatment can be carried out without these simplifying 
assumptions; however, the final expressions will be much more cumbersome 
and the number of polarizability terms characterizing a bond will be larger. 
In the following we will restrict ourselves to non-cyclic molecules with bonds 
having cylindrical symmetry. 

CALCULATION OF dau,/dQ* 

Let us consider a non-cyclic molecule with N bonds. In the normal coordinate 

treatment internal coordinates are used (bond stretchings, interbond bendings, 


etc.), because in terms of them the potential constants have a clearer physical 
meaning. Then, symmetry coordinates are formed as linear combinations of 
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internal coordinates, which belong to the different irreducible representations 
of the point group of the molecule; thus the secular equation is factored out 
into several equations of lower order. However, a tensor referred to space- 
fixed axes (fixed as far as the vibration is concerned) such as the polarizability 
tensor cannot be defined in terms of internal coordinates. Wolkenstein and 
Eliashevich, and Long make use of some sort of intermediate coordinates. 

Long defines the so-called intensity coordinates. For the molecule of NV 
bonds which we are considering, there are N bond stretchings, and 3N direc- 
tion cosines of which only 2N are independent because of the normalization 
condition. For the component a,,, the N bond stretchings and the 2N direction 
cosines cos(nu), cos(nv) are used as intensity coordinates. 

The following nomenclature will be used: 


R‘* = ith internal coordinate, 
S’ = rth symmetry coordinate, 
Q* = kth normal coordinate, 
J? = pth intensity coordinate, 
d, = displacement vector of the ath atom during vibration. 
The symmetry coordinates are related to the normal coordinates by means 
of the following matrix defined by Wilson 


[8] ST = x LiQ* 
and 
[9] : Q = dX ares, 


where LZ and L~ are obviously reciprocal matrices. The use of superscripts 
and subscripts has been kept here for consistency. 
The symmetry coordinates are defined from the internal coordinates by 


means of 
[10] Sa Fur. 
t 


If there are redundancies, i.e. dependence between the internal coordinates, 
the matrix U will be rectangular; however, by making use of the redundancy 
relations it is possible to solve for the internal coordinates in terms of the 
symmetry coordinates, so that 


[11] R‘ = > Tis". 


Since the symmetry coordinates are orthonormal, including the redundancies, 
the matrix T is the transpose of the matrix U. 

For the part involving the direction cosines, the displacement vectors d, 
can be used as intermediate coordinates. The displacement vectors cannot be 
written directly unless explicit use is made of the Eckart conditions (3), since 
there are six more displacements than internal coordinates. On the other hand, 
the internal coordinates can be written as linear combinations of the displace- 
ment vectors as 


[12] R' = > sj-d,, 
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where 8; is the s-vector associated with the ith internal coordinate and the 
ath atom as introduced by Wilson. Explicit formulas for the 8-vectors corre- 
sponding to the different types of internal coordinates have been given in the 
literature (2, 4, 5, 9). By taking into account the consistency of Wilson’s method 
with the Eckart conditions by which the latter are automatically taken 
care of (7) the displacement vectors can be expressed as 


[13] ds = wa Dy 82gisR' 
F 


where u, is the reciprocal of the mass of the ath atom and g,; is the kinetic 
energy matrix defined by 
27 = > giyR'R’. 
ij 

Eq. [13] will be derived in the Appendix and corresponds to an expression 
given by Crawford and Fletcher (1) and used subsequently by Long (6). 

We can now obtain the elements of the derived polarizability tensor, 
namely, 


OO _ a Oetur al’ 


co al’ ag 
Now, WN of the intensity coordinates are bond stretchings which may be 
designated by the indices p = 1, 2,..., N; the rest designated by the indices 
bp = N+1,...,38N are direction cosines. Then 
Datyy SR Dyn OI? | wR Day, OI” 
14 eas Seas 2S = =. 
14] aQ* d, al? aQ* hee al’ aQ* 


The first term on the right-hand side of Eq. [14] can be obtained very easily. 
From Eggs. [8] and [11] we have 


[15] Ri = Soringe’. 
rk 
and if we designate the internal coordinates corresponding to bond stretching 
by q", then 
g" = 2) TrLi". 
rk 
Then from Eq. [7] 


Yi Betas BT” _ 5 Sawn dg" 
[16] = oF 3 dq” ad 
Oein2 


0 n Tr 
ee » | cos(nu)cos(n 5 (ou1 ans) [recs 


Although the index r in Eq. [16] formally runs over all the symmetry coordi- 
nates, 7} is clearly zero unless r refers to a symmetry coordinate involving bond 
stretchings. 

For the second part on the right-hand side of Eq. [14] we have: 


3N Boru» ol” > Doe ol” dd, oR! 


[17] p=N+1 ol? a0" . p=N+1 al’ ai dd, AR’ ag 
3N 0a P . i 

uD : 

= > —); ad, Ste8ssT Ln, 


p=N+1 0 ajri 
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where use has been made of Eqs. [13] and [15]. The index 7 runs formally over 


all internal coordinates. 
We can define a set of vectors Sj each of which is associated with an atom 


and a symmetry coordinate by means of 


[18] e= DTS 


in terms of which the displacement vectors can be written as 


[19] d, > o> SiG,.S* 


which is similar to Eq. [13] and where G,, is the kinetic energy matrix defined 
by 


[20] oT = > G,,5S'S', 
rt 


which is related to g;; by means of 
[21] Gre = Do gusT iT? 
ty 
From Eq. [18] we have 
ol” j mrirr ol” tori irr 
D ad, Sabai! Li — Zz ad, 38! Mek isT Le 
ar’ 
atr dd, 
al’ - 
zz ad, Suta(L ‘Ys 


at 


[22] _ . SiuaG erLt 


where the last step is obtained from the relations between the matrices L 
and L~! introduced by Wilson (9). 
From Eqs. [7] and [22], Eq. [17] can be written as 


3N n 

Day, OL 

9 7 GGus. CF 
P| x al? aQ* 


= D> (an1—an2) » | on(n0)<2-cos(nu)-+cos(n)-2-cos(ms)-Shue(Z-% 


In the summation with respect to the index a involved in Eq. [23] only two 
terms will appear, those corresponding to the two atoms which form the bond. 
If we call these atoms 7 and j, then 
0 . 0 . I 
ad cos") = —sin (nu) 5g (nu) = sin(nu)fixQn, 

where f?, is a unit vector perpendicular to the ji-bond towards the u-axis, 
and Q, is the reciprocal of the bond distance. This result follows immediately 
from the meaning of the partial derivative as the gradient of the angle at 
the equilibrium position (see Fig. 1). The vector f/, can be written as a linear 
combination of a unit vector from the jth atom to the ith atom e,,, and a 
unit vector along the u-axis e,, 
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Fic. 1. The gradient of the angle between the ji-bond and the u-axis. 


fi, = Satna lx@11008(00)1 
Therefore 
[24] or = -—» copies) = Q,[e,—e,cos(nu)]. 
ad, ad, 
The final expression for the derived polarizability is then 
[25] a = a [cos nn)eos ne) 2 (a —ons)+0u-20 raz 


= 7. (an1—On2) Q,[cos(nv)e,+cos(nu)e,— 2cos(nu)cos(nv)e 51] 
-(uSi—wjSj) (Li. 


From Eq. [7] it is clear that in the trace the second part of Eq. [25] contributes 


from the normal property of direction cosines. Therefore only the first term of 
Eq. [25] contributes to the trace and from Eq. [16] 
[26] ie = Zs Fo (cier+ Pees) THEE. 

w OQ “ar 8g 
For a symmetry coordinate of totally symmetric type, all the elements 7” 
are the same, since symmetry coordinates are formed as linear combinations 
of equivalent internal coordinates; by using now the orthogonality property 
of symmetry coordinates it is seen that the trace is zero for any vibration 
other than the totally symmetric ones and therefore the corresponding 
depolarization factor is 6/7. This has been shown similarly by Wolkenstein 
and Eliashevich (10) and by Long (6). 


AN EXAMPLE 


As an example we will consider the tetrahedral CX;Y molecule, shown in 
Fig. 2. The direction cosines are given in Table I. The S vectors can be ob- 
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The tetrahedral CX;Y molecule. 


Fic. 2. 
TABLE I 
DIRECTION COSINES FOR THE TETRAHEDRAL CX3Y MOLECULE 








Axes 
Bond ee = are area semana 
x y z 
cy 0 0 1 
CX; 0 sin @ —cos 6 
ex, —}4/3 sin 0 —}sin@ —cos 0 
CX; 34/3 sin 0 —} sin 0 —cos 0 





tained from Eq. [18] and the s vectors given by Wilson. They are given in 
Table II, where the unit vectors are shown in Fig. 2 and the Q’s are the 
reciprocals of the bond distances in the equilibrium position. 


TABLE II 
S VECTORS FOR THE TETRAHEDRAL CX3Y MOLECULE 


Symmetry coordinate 





Atom —_—————e 
1 2 3 
7/: ‘ a. 

Cc —=@ — sa lertertes) [ $70.- a, | (e:te:+e,)— V (80:+Q,)e 
Xi 0 Fi e; = ve Qs(e:-+3e2-+3e; —3e) 

Xe 0 a e» _ “ Q,(e2+3e; +3e; —3e) 

X; 0 a es = 4 0. (es +3e: +3e: —3e) 

¥ e 0 4 Q,(e: +e2+e3+e) 
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From Eggs. [25] and [26] and use of the values given in Tables I and II, 
we obtain 


Oar, , 


ag = aioli t+3 ff (a1 — az) +3e2)Li 
+3(3-*) (az1—arge) Oz {2[3-* (LZ) §—3(L")$+-302(L") $c tuxQ2(L*)§}. 


The trace and the anisotropy are then 


Ay = (aj: +2a40) Li +3? (ab1+20%2) Li 


and 
By = (a1 —e}2)L4—3 (ah — a2) Li 
—4 (3) (@s1—az2) Qz{2[3 (LZ )i-3(L") 3 +302(L" suc tux Or(L™ 5}. 
The prime indicates the derivative of the bond polarizability with respect to 
bond stretching. 
APPENDIX 


In a previous paper (7) the matrices D and Q were defined as relating the 
internal coordinates to the mass reduced displacement coordinates. By Eq. 
[10] of Reference 7 we have 


> GDit+ dX GD? = &. 
j Dp 


Multiplying by Di 6° we get 
2 GDDs"+ 2) QDDs” = 2) DS 
cb 


jbc poe 
from which, by use of Eqs. [12] and [14] of Reference 7, we obtain 


> Ce! = D;. 
J 


Multiplying now by gi and summing over 7 we finally obtain 


i 
Or = i Dag ix 
t 
By splitting the index a, the last equation can be written as 


9 1 i 
“aRE = ar Digi = ae Sn Og ix- 


By combining now the three components of d,, the displacement vector of the 
nth atom, we get 


d, oa Hn De SZ nR", 


which is Eq. [13] of this paper. 
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Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


Atomic Masses of Ru®, Ru, Os!8°, and Os!* * 


The atomic masses of Ru, Ru®, Os!8%, and Os!* have been measured using the Dempster- 
type mass spectrograph at McMaster University. Both the instrument (2) and the method (3) 
have been described in detail previously. The mass differences obtained in this work and the 
atomic masses deduced therefrom are given in Table I. 


TABLE I 
MEASURED MASS DIFFERENCES AND COMPUTED ATOMIC MASSES 





Doublet Mass difference (m.m.u.) Atomic mass (a.m.u.) 
403189 — Cus 56.6 +3 Os!89 = 189.0177 + 9 
4Os!92 — Nis 59.9 +2 Os!® = 192.0224 + 8 
sPt1%s — Ru% 770+3 Ru% = 97.9363 + 4 


98.9378 + 4 


$Pti98 — Ru? 76.3 + 2 Ru 





In order to calculate each of the above masses it was, of course, necessary to know the mass 
of the other member of the doublet. The values of these reference masses were taken to be: 
Cu® = 62.941926 +6 (1), Ni® = 63.94755 +7 (4), Pt!%* = 196.0267 +6 (5), and Pt! 
= 198.0283 + 6 (5). In the case of Pt!%8, an error in arithmetic leading to the wrong value 
given in Reference 5 has been corrected. 
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CORRECTION 


Vol. 30, 1952, pp. 574-575. The measured values of the electric field in Figs. 2 and 3 are in 
error and should be multiplied by 0.68 to obtain the correct r.m.s. fields. 
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